NASA CR-132528 


APPLICATION OF THE SCALAR AND 

i 

VECTOR POTENTIALS TO THE 
AERODYNAMICS OF JETS 


by 


H. L. Russell and J. G. Sklfstad 


\ 

August 1973 


t ; Thermal Science and Propulsion Center 

School of Mechanical Engineering 
Purdue University 
Lafayette, Indiana 


NASA Grant No. NGL 15-005-094 
Purdue Research Foundation PRF-5551 


(NASA-CR-13252 8) APPLICATION OP THE N74-34702 

SCALAR AND VECTOR POTENTIALS TO THE 
AERODYNAMICS 0? JETS (Purdue Univ.) 

142 p HC $10.25 CSCL 2CD 



w 


Unclas 

G3/12 52010 



ACKNOWLEDGMENTS 


This investigation comprises part of a research program on 
the aerodynamics of jets, sponsored under NASA Grant No. NGL-15- 
005-094. The authors gratefully acknowledge the support and 
especially the encouragement and cooperation of the project monitor, 
Mr. Richard J. Margason of the NASA Langley Research Center. Dis- 
cussions with Dr. R. E. Lynch, Div. of Mathematical Sciences, 

Purdue University, regarding numerical methods were most helpful. 

The work was completed in 1971; publication of this report has 
been delayed largely because of other commitments of the senior 
author (J.G.S.). 



in 


41 


preceding page 


blank not filmed 


TABLE OF CONTENTS 


LIST OF TABLES . . . 
LIST OF FIGURES . . 
LIST OF NOMENCLATURE 

abstract 

I. INTRODUCTION . . 

II. ANALYSIS . . . 


Page 

v 

.vii 
vlii 
• xii 
. 1 
. 3 


III. 


IV. 


A. Vector Field Theory 3 

1. Fundamental Theorem 3 

2. The Scalar and Vector Potentials !!!!!!!!!! 4 

3. Field Components on the Boundary .... in 

B. Solenoidal Fields ;V 

C. Time Dependent Fields !!!!!!!!!!!! ‘ io 

D. Axisymmetric Fields 13 

E. Green's Function for a Semi -infinite Domain ....... 13 

THE AERODYNAMIC FIELD OF A JET 20 

A. Method Based on the Scalar Potential. ... 22 

B. Method Based on the Stokes Potentials ! ! ! 24 

C. Vorticity in the Jet Flow . .... 25 

THE AXISYMMETRIC JET IN A QUIESCENT ATMOSPHERE 37 


A. Formulation of the Problem 

B. Numerical Method !.!!!!!!! 

C. Computations for Selected Vorticity Distributions ! ! ! . 

1. Vorticity Distributions Corresponding to 

the Selected Velocity Profiles 

a. Mixing Layer in the Potential Core .....!] 

b. The Developed Region 

2 . Regional Contributions to the Aerodynamic Field ! ] ! 

a. Relative Contributions of Different 

Regions of Vorticity along the Length 
of the Jet 

b. Contributions from Radial Elements of the 
Distribution of Vorticity 


41 

47 

55 

56 
56 
60 
67 


67 


i 


• 67 



3. The Potential Field 70 

4. Summary of Results for the Selected 

Velocity Profiles 73 

D. Approximate Computational Models 76 

1. Triangular Vorticlty Profile 81 

2. Uniform Vorticlty Profile 83 

3. Vortex Sheet Model 83 

4. Results of Computations Based on the 

Approximate Models 84 

V. DISCUSSION 89 

A. Contribution of E Relative to U<. 89 

B. Comparison of Computed Results T 

with Experimental Data and Other Theories 90 

1. Comparison with Experimental Data 90 

2. Comparison with Other Theoretical Models 94 

C. Comparison of Computations for the Approximate 
Models with the Computations for the Reference 

Model (Squire and Trouncer Velocity Profile) 97 

D. Evaluation of the Integral Moment of Vorticity 99 

VI. CONCLUSIONS 105 

LIST OF REFERENCES 108 

APPENDICES 

Appendix A: Annotated Bibliography on Axisymmetric Jets... 110 

Appendix B: Derivation of the Equations for the Aerodynamic 

Field of an Axisynmetric Jet 123 

Appendix C: Algorithms for Complete Elliptic Integral 

Functions 129 



V 


v 1 

LIST OF TABLES 

Table Page 

1. Field Point Velocities - Contributions from 
the Mixing Layer of the Potential Core 
[Squire and Trouncer (cosine) Velocity 

Distribution] 58 

2. Field Point Velocities - Contributions from 
the Mixing Layer of the Potential Core 

(Abramovich Velocity Distribution) 59 

3. Field Point Velocities - Contributions from 
the Developed Region [Squire and Trouncer 

(cosine) Velocity Distribution] 64 

4. Field Point Velocities - Contributions from 
the Developed Region (Gaussian Velocity 

Distribution) 65 

5. Field Point Velocities - Contributions from 
the Developed Region (Schlichting Velocity 

Distribution) 66 

6. Field Point Velocities - Effects of Extending 
the Region of Integration [Squire and Trouncer 

r (cosine) Velocity Distribution] 68 

7. Field Point Velocities - Contributions from 

Radial Elements of the Jet 69 

8. Field Point Velocities - Contributions from 

the Potential Field 71 

9. Total Field Point Velocities [Squire and 

Trouncer (cosine) Velocity Distribution] 74 

10. Total Field Point Velocities (Gaussian 

Velocity Distribution) 77 

11. Total Field Point Velocities (Schlichting 

Velocity Distribution) 79 

12. Field Point Velocities - Contributions from 

the Combined Mixing and Developed Region 87 

I 



vi 

7 a Lit Page 

13. I xporii itntal and Theoretical Values of -C x 10 A 91 

P 

14. Comparison with Other Theoretical Models - C x 10^ J(> 


LIST OF FIGURES 


vi i 

r,flur? Page 

1. Geometry Associated with the Vector Notation . 6 

?. Cylindrical Coordinates 

3. Spherical Coordinates 

4. Semi-Infinite Domain ig 

5. Aerodynamic Field and Jet Flowfield 21 

6. Domains for Aerodynamic Field 23 

7. Peripheral Component of Vorticity 27 

8. Small Section of Jet Flow 29 

9. Geometry for the Line Integral 

10. Axisymmetrlc Turbulent Jet 38 

-► 

11. Geometry Associated with the Volume Integrals for U c . . . 44 

bT 

12. Geometry Associated with the Surface Integral for ? 46 

13. Geometry of the Jet Flow Domain and Transformed 

Jet Flow Domain 

14. Velocity and Vorticity Distributions in the 

Developed Region 52 

15. Approximate Models for the Vorticity Distribution 85 

It). Comparison of Experimental Data with Theoretical Values. . 


92 


vm 


LIST OF NOMENCLATURE 

A exact value of an integral 

-> 

A vector potential 

A-j surface area 

- ► 

A] a vector normal to and equal in magnitude to the surface area Ai 

C a constant in the expression for velocity potential. Equation (79) 

C p pressure coefficient 

e virtual origin of the jet 

E (k) complete elliptic integral of the second kind 
E, (l mass entrainment rate 

-> ^ *>• 

Ep irrotational-solenoidal contribution to the vector U from the 

surface integrals Equation (13) 

0 Green's function of the second kind 

II a harmonic function 

1 a numerical approximation for an Integral 

I m integral moment of vorticity 

k modulus for complete elliptic integral functions 

K-j a constant in the expression for the mass flow rate. Equation (39) 

K (k) complete elliptic integral of the first kind 

M momentum flux 

m exponent in the empirical expression for the velocity distribution 

in the jet 

iii dipole strength 

n the unit outward normal vector of a surface element 

p pressure 



lx 


> 



p. 

*> 

Q 


R 

R 


r 


«V 

2 

s 

dS 

d$ 

T 

-+ 

U 

u 

m 

-► 

W 

->• 

X 

X, 

x, 

z 


r 


position vector for an image of a source point 

ambient pressure of the atmosphere 

a vector field such that Q = u - T 

distance from a source point to a field point 

magnitude of the distance between a source point and a field 
point 

radial coordinate in a cylindrical coordinate system 
inner boundary of the mixing layer of the jet 
outer boundary of the jet 

radius at which the velocity is half that at the axis of the jet 
a surface bounding a region in space 
a differential surface element 

the product of the unit outward normal vector and the differential 
surface element 

a vector field such that vxT=vxU,v-T=v*U, and 

-*•-*-*•-*■ 

T * d$ = U * dS 

a vector field 
maximum axial velocity 
a vector field such that v 2 W = - 0 
a position vector 
Y, Z rectangular coordinates 
y, z transformed coordinates 
axial coordinate in a cylindrical coordinate system 

Greek Symbols 


6 


circulation 
the delta function 

the error of a numerical approximation 


f i a transformation variable 

c: a cylindrical coordinate 

o, p, 4> spherical coordinates 
a transformation variable 
y a transformation variable 

f (18.5 r-j )/z-j for the axlsymmetric jet 

n mass density 

t a region in space 

4. scalar potential 

il> gradient of a scalar 

-► 

vorticity vector field 

Subscripts 

C pertaining to the centerline of the jet 
I Irrotatlonal field 

j pertaining to the jet 

m maximum value 

p peripheral component 

r radial component (cylindrical coordinates) 

S solenoidal field 

x transformed radial component 

y transformed axial component 

z axial component (cylindrical coordinates) 

r peripheral component (cylindrical coordinates) 

c azimuthal component (spherical coordinates) 

P radial component (spherical coordinates) 



xi 


i a region in space 

'!' colatitude component (spherical coordinates) 

1 pertaining to a source point 

Superscripts 

C core region of the jet 
D developed region of the jet 

" (Prime) per unit length 



xn 


ABSTRACT 


An analytical investigation was conducted to examine the applica- 
bility of a method based on the Stokes potentials (vector and scalar 
potentials) to computations associated with the aerodynamics of jets. 

The principal merits of the method were found to be that the aero- 
dynamic field near the nozzle could be well represented and that 
the influence of a nonuniform velocity profile at the nozzle exit 
plane could be explicitly determined. Other aspects of the computa- 
tions indicated accuracy comparable with other methods but at the 
expense of more complexity in the computations. An additional 
benefit of the theory was that it provided the rationale for de- 
veloping useful approximate analytical models for computations of 
the type considered. 

Computations were made for an axisymmetric jet exhausting into 
a quiescent atmosphere for the purposes of exploring the computational 
aspects of the method and for comparison with the available experimental 
data and other theories. It was found that the velocity profile at 
the exit of the nozzle and the integral moment of vorticity in the 
jet flow were the most significant factors in the computations. The 
latter is shown to be in agreement with the asymptotic form of the 
theory. Those factors suggest that knowledge of the velocity at the 


x 1 i 1 


* x ' s of the jet, together with the location of the half-velocity points 
along the jot are sufficient to yield reasonably accurate aerodynamic 
field computations for the particular flowfield considered. Computa- 
tions support that conclusion. 

One of the approximate analytical models considered in the study 
essentially replaced the jet flowfield by a conical vortex sheet. 

That model was the simplest considered and yielded accurate results. 

It is, accordingly, recommended for computational purposes. 

The comparison with experimental data showed agreement within 
the scatter of the available data. The theory indicated clearly the 
inconsistency among the different theoretical characterizations of 
jet flowfields, and suggests the need for more accurate relations. 


I. INTRODUCTION 


A number of problems in fluid mechanics involve the interaction 
between viscous and inviscid flows. For most purposes, as in boundary 
layer theory, the inviscid flow is determined without consideration 
of the viscous flow region whereas the effects of the Inviscid flow 
on the flow in the viscous layer are considered. More recent interest 
in the aerodynamics of jets, the "strong interaction" problems, separat- 
ing and reattaching flows, etc., have placed more emphasis on the 
mutual interactions of the two regions. Analytical models for such 
flows are in various stages of development. One analytical method, 
applicable to incompressible flows or for aerodynamic purposes, 
based on the Stokes potentials (scalar and vector potentials), has not 
yet received much attention. 

It is the purpose of this Investigation to examine the latter 
method in the context of problems associated with the aerodynamics of 
jets. Such flows are of considerable significance in the aerodynamics 
of V/STOL aircraft, the dispersion of pollutants, etc. Analytical 
methods for handling them are not yet well developed. The emphasis 
in this connection is on the influence of a turbulent jet in inducing 
an aerodynamic flow. The flow in the turbulent jet is presumed known 
either from experimental data or from other theoretical considerations. 

As part of the investigation, the properties of the turbulent flow 
most influential in computing the aerodynamic field were evaluated. 
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Asymptotic relations were developed to yield greater insight into the 
method. In the interest of simplifying such computations, several 
approximate models were examined. Numerical computations were made 
for the case of an axi symmetric jet exhausting through a plane baffle 
plate, serving to further classify the advantages and disadvantages of 
the method. 

The use of the scalar and vector potentials arises naturally from 
the fundamental theorem of vector field theory. Most applications 
have been associated with electromagnetic field theory; applications 
in fluid mechanics have been rare. Since the concept is not familiar 
to those with prime interests in fluid mechanics, a brief review of 
the relevant aspects of vector field theory have been included In 
the text. 

The results of the investigation show that the method offers 
unique advantages in the computation of the Induced aerodynamic field, 
particularly close to the origin of the jet. At more distant posi- 
tions, the method yields results comparable to those obtained by other 
less complicated methods. An additional benefit provided in the appli- 
cation of the method is that It offers an Independent means for com- 
paring the consistency of different representations of the turbulent 
jet flow. The examination of approximate models for the jet flow 
showed that quite reasonable approximations for the induced aerodynamic 
field could be obtained from knowledge of only the centerline velocity 
and the positions of points in the jet at which the velocity is half 
that on the axis. 



II. ANALYSIS 


Since the basis for the analytical methods to be considered lies 
in vector field theory, an Introductory description of the relevant 
elements of that theory is presented in the first part of this section 
The remaining parts of this section present discussions of solenoidal 
fields, time dependent fields, axisymmetric fields, and Green's 
function, respectively. 

A. Vector Field Theory 
1. Fundamental Theorem 

“► 4 -*• 

Consider a vector field U (X) which depends upon the position X 

in space. The field U may be uniquely defined in a simply connected 

region t by specifying the divergence and curl of U throughout t 

and the normal component of U on the surface S bounding the region t. 

1 “► 

A proof of uniqueness, following Lass , is as follows: Let T be 

another vector field with divergence and curl equal to those of U, 

and with T * dS = U * dS on S.* Select a vector Q which is the 

vector difference between U and T, Q = U - T. This vector is 

*> 

irrotational , solenoidal, and its scalar product with dS is zero. 


*dS = n dS, where n is the unit outward normal vector of the surface 
element, and dS is the differential surface element. 



The vector Q, since it is irrotatlonal , may be represented In terms of 
a scalar function <t> such that 

<$ = - v <j> 

The solenoidal condition then yields 

V*Q = -v 2 * <f, = 0 


Applying Green's theorem in the first form yields 


J 


I 


9 0 

[41V <f> + (v<f>)^] d T = I + v« * dS - I 4. Q • dS * 0 


L 


and therefore the volume integral of ( v<t> ) is zero. Hence, V 4 > = 0 

-> ■> -> 

inside t, and Q =-V 4 > = 0, so that U = T Inside t. For an Infinite domai 
the condition that U ->0 as 1/R* + G as R -*■ • with e > 0 is required. 2 


2. The Scalar and Vector Potentials 
Explicit representation of the vector field in terms of the volume 
distribution of its divergence and curl, and its normal component on 

3 

the surface may be obtained by Introducing the Stokes potentials , 
normally called the scalar and vector potential of the field. Follow- 
ing Morse and Feshbach 4 ^, let 


W (X) = 


1 

~ET 



u 1 (x t ) 

R (x, Xj) 



( 1 ) 
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where R = |X - x | 

X * the field point position vector 
Xj * the source point position vector 
dxj ■ the volume element with respect to 
integration over the source field 

Refer to Fig. 1 for the geometry associated with this representation. 
Equation (1) represents a solution of the vector Poisson equation 
v 2 W = - U. Thus, 

U = - v 2 W = - v (v • W) + 7 x (v x W) ( 2 ) 

Equation (2) yields a representation of the vector field U in terms 
of the negative gradient of a scalar and the curl of a vector. Let 


and 


<M V • W, 

-► “► 
A = 7 x W 


(3) 

(4) 


The scalar <p and the vector A are the scalar and vector potentials of 
the vector field U, respectively. The general representation of a 
vector field as the sum of an irrotational field and a solenoidal 
field has been attributed to Helmholtz. Denoting the irrotational 
field by Uj and the solenoidal field by 11$, then, 

1 * J I + “s 

where Uj * - 7 <p 

U<- = 7 x A 
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FIELD POINT 



FIGURE I. GEOMETRY ASSOCIATED WITH 
THE VECTOR NOTATION 
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Taking the divergence of Equation (1) with respect to the field 
► 

point coordinates X yields* 



dt 


1 


where the operator v is with respect to the field point X. Because of 
the symmetry of R with respect to X and Xj, however, one has 


7 <p> - - 7 i <if> 


( 5 ) 


where the operator is with respect to the source point X^. Thus, 



1 1 I “► 

Noting the Uj • Vj (^) = v i * (pr U ^ . Uj and making use of the 
Gauss divergence theorem, the scalar potential may be finally written 


-+ * -► 



Equation (6) relates <t> to the divergence of U in t and to the normal 

~y 

component of U on S. 


♦Variables with a subscript 1 denote evaluation at the source point Xj. 
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Taking the curl of Equation (1) yields 


kl ^1 x 7 W dx l 


Using the identity. Equation (5), results in the expression 


h j x U 1 x 7 1 $ dT l 


By employing the vector identity 


5 1 x 7 1 ( K> * R 7 1 x U 1 


Vj x ( r ). 


the integrand may be expressed in the form 


A = - 


4* J T r 1 * 'R 


U, V, x U, 
Vi x («i) - -hr-1 


Making use of the appropriate vector integral theorem related to the 
Gauss theorem, one has 


V 1 x ( R~ } dT l = 


Uj x dSj 


which yields, finally 


7 1 X U 1 ri ,1 f U 1 X dS l 

“I! — dt l + I — R 


Equation (7) relates the vector potential to the curl of U throughout 
■( and to the tangential components of U on S. 

Taking the negative gradient of Equation (6) with respect to the 
field point X, and using the identity 


one obtains 


. v ^ ^ ( V] • U,) 


J dT l 


h (U 1 


n ^ R 
Hi ) ~*r 


3 dS l 


( 8 ) 


Equation (C) shows the irrotational field is related to the divergence 
of U throughout t and to the normal components of U on S. 

Taking the curl of Equation (7) yields 


v x A = 


1_ 
4 71 


(v-i 


x U-j ) x 


-> 

R 


3 dl l 


+ k < U 1 


x n.j ) x 


R 


3 dS l 


(9) 


Equation (9) shows the solenoidal field is related to the curl of U 
throughout t and to the tangential components of U on S. 

Equations (6) and (7), and (8) and (9), respectively express 
the potentials and the field values In terms of the volume distribution 
of the divergence and curl of U and the components of U on the boundary 
of the domain. 
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3. Field Components on the Boundary 
The representation given by Equations (8) and (9) does not 
correspond with the statements made In Sec. II A 1. In the sense 
that the surface integral in Equation (9) Involves the tangential 

-V ->■ 

components of U on S rather than the normal components alone. It 
may be shown, however, that the combination of the surface integrals 
in Equations (8) and (9) give rise to a field that Is irrotatlonal 
and solenoidal, so that specification of the divergence and curl of 
the field throughout the region, together with only the normal compo- 
nent of U on the surface of a simply connected region Is sufficient 
to define the field to within a constant vector. 

To demonstrate the latter, consider the vector field to be repre- 
sented by the sum of three velocity fields in the region t and on Its 
surface. 

3 ■ 3 i,* K * 3 p (10) 



Let the region exterior to t and extending to Infinity contain no sources 
for the field, that Is, v • U 8 v x U * 0 outside t. Since the field 



11 


U in t due to a given distribution of sources in space does not depend 
on the particular choice of x, let a new domain x"be defined to include 
t and part of the region exterior to x. For any point in x, then, 
the volume integrals yielding U Jt and U St do not depend on the choice 
for t\ Thus, the remaining field E p and therefore the surface 

-f 

integrals which determine E p must be independent of the choice of t r 
Taking x' sufficiently large, the surface integrals vanish subject to 
the condition that U -*• 0 sufficiently fast 2 as R -* and so they 
vanish identically. Thus, E p = 0 for this case. For the general 

situation, then, the surface integrals may be seen to represent the 

*> -► -+ 

contribution to U of all sources (v * U, v x U) outside x. Invert- 
ing the domains in the foregoing example, let v*U=vxU=0 
within x, and non-zero in some bounded region exterior to x. Thus, 

E p is then solenoidal and irrotational within x and may be represented 
by a harmonic function given in terms of its normal components over S, 

tp = - L v G V U, * dS, (14) 

*> -f L 

where G (X, ) is Green’s function of the second kind for x. This 

problem of finding a potential function, harmonic in a region x, 
bounded by a surface S, and having normal derivatives given on the 
boundary is known as Neumann's problem. The general representation 
of the field as given by Equation (10) is also described by Batchelor^ 3 ) . 

B. Solenoidal Fie lds 

Most of the discussion in the following section is concerned with 
incompressible flow, for which the divergence of the velocity field is 
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zero. Ihus, the volume integral involving the di vorgence of the 
souroe Mold Iquation (II) for this case will he iv.ru. Ihe velocity 
field may then be expressed as 


u * u s, 



00 


Thus, the determination of a solenoidal field is completely determined 
by evaluating a volume integral of the curl of the vector field and 
then solving the appropriate problem in potential theory. It is this 
problem that is of principal concern in the following sections. 


C_. Time Depe ndent Fi elds 

In many situations in fluid mechanics, time-varying phenomena 
must be considered. Frequently, this phenomena is in the form of 
turbulence. As long as one is not concerned with relativistic theory 
and relativistic effects, time may be regarded merely as a parameter 
of the vector field. For that case, all of the above relations are 
valid for a time varying field U (X, t). Because of the linear char- 
acter of the vector equations, the time averaged field is simply that 
computed from the time averaged divergence, curl, and surface normal 
components of the field. Stating this another way, the time averaged 
field is that field computed from the divergence and curl of the time 
averaged field and the surface normal components. Thus, time averaged 
equations have the same form as the steady state equations. 
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l 1 . Axi symmetric Fields 

For axisymmetric fields with no tangential component 




(cylindrical coordinates) 


or 


then 


U = U i + U i 

p o <t> 


- 3 U 3 II 

v x U = ( - r z 


(spherical coordinates) 


3 z 3 r ' o 


- ) i n (cylindrical coordinates) 


or 


and 


1 ^ P U 3 U + 

y x U = — ( — — - 1 — ^-) i 0 (spherical coordinates) 


S T r + x o T e + X z V - t x rl T rl * x el T el + X zl T zd 

(cylindrical coordinates) 

R = [X T + X T + X il - [X , T , + X,, T . + X , T q1 ] 

r> r> 4> <J> o o J pi pi <\>\ 1 ol ol J 


(spherical coordinates) 


The geometry of the fields are shown in Figs. 2 and 3. The coordinates 
with subscript 1 correspond to a source point and the nonsubscripted 
coordinates correspond to a field point. 


E. Green's Function for a Semi -Infinite Domain 
Assuming the divergence, the curl, and the normal component of a 

-V 

vector field U are specified within a given domain, it can be useful 

to employ Green's function of the second kind to determine the 

■> 

component of the vector field denoted by Ep (Equation 14). Since 
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FIGURE 3. SPHERICAL COORDINATES 
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E p is irrotational it may be represented by the gradient of a scalar <|*. 
Thus, it follows from Equation (14) that 


♦ 




G ( X , X-j ) U-j 





2 

with v i|> = 0 throughout the domain. 

Consider a simply connected finite region t. To obtain the Green's 
function for this region, the following vector function is used. 


*1 v iT‘ (T v *1 


where ^ is harmonic within the domain t. Express Green's theorem in 
the second form 


I 


(*i v' 


l 

R 


1 

IT 


^ ) 


dx, = 


i 


u 


1 an 




1 

IT 


3*! 

arT 


-) dS 


l* 


(16) 


The volume integral involving v <|^ is zero. Solving the left side 
of Equation (16) for field points within the region t, it follows that 

<’ 7 > 

•'S 

o 

since v 1/R = 0 except at R = 0 where it takes on a value of - 4n. 
Equation (17) appears to indicate, that to obtain a solution for t|», 
both i|>.| and a^/an must be known on the surface. This is contradictory 
to Equation (14) which requires only the normal derivative (U^ • dS^ 

aiK 

= — - dS^). It is possible to express the integrand of Equation 07) 
in terms of a^/an alone. 
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In accordance with the presentation by Kellogg ', consider a func- 
tion II which is harmonic in x, and form the vector function ^ vH 
- H analogous to that shown above for H = 1/R. Thus, Green's 
theorem in the second form for this vector function yields 



Addition of Equations (17) and (18) yields 

*--lf [ *i lr ( i + H) ■ sir H)J dS r («) 

The combined scalar (^ + H) is called Green's function of the 
second kind. The scalar potential may be eliminated from the 
integrand by requiring that the normal derivative of Green's function 
on the surface be zero. 

The function H must differ from - 1/R, in general since H is 

harmonic at all points within t whereas 1/R has a singularity within x. 

a 1 

It may be noted that the surface integral of — is -4tt and the 
surface integral of the normal derivative of the harmonic function H 
is zero. If 3H/3n differs from - ^ r b y a constant at points on the 
surface, the surface integral involving yields a constant. With 
that choice for aH/an , Equation (19) provides a relation for <p in 
terms of a^^/an on the surface to within a constant. The determination 
of Green's function G (X, X 1 ) then reduces to solving the following 
problem. 

G (X, = 1+ H (X, X,) (20) 

ORIGINAL PA-GEE 

OF POOR Q UAU1 
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where H is a hannonic function in x satisfying the boundary conditions 


ill 

an 


3 1/R (X. X,) 
an 


+ C 


( 21 ) 


on the surface of the domain where C is a constant. 

As an example, consider the semi-infinite domain shown in Fig. 4. 
The point P 1 is the image with respect to the plane Z = 0 of the 
source point r^. If H is given the harmonic value 1 / | r - | , and 

substituted into Equation (21), the results obtained are 


Thus 


3 1 




on S, 


( 22 ) 


G = — L— - + 7--V- (23) 

I r - r ] | |r - P ] | 

-* -* -+ 

Note that the quantities |r - r-j | and |r - P^| are equal whenever r^ 
is on Z = 0. The required Green's function for the semi -infinite 
domain is then simply 2/R. For an infinite domain the appropriate 
Green's function degenerates to 1/R. Green's functions for various 
domains are discussed at length in Reference 4 (b). 
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FIGURE 4. SEMI-INFINITE DOMAIN 
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III. THE AERODYNAMIC FIELD OF A JET 

Referring to Fig. 5, a jet exhausting Into an airstream has the 
property that most of the turbulent flow remains confined within a 
relatively narrow, curved column. There generally will be a turbulent 
wake behind the jet, but for this discussion. It Is convenient to 
consider the entire flow field to be comprised of two parts: the 

i±t flow , and the flow outside the jet, or the ae rodynami c field. 

The jet flow occupies that region containing any significant vorticity. 
The vorticity in the aerodynamic field, then, is presumed to be neg- 
ligible. Experimental data for jet flows support such a separation 
of the flowfield based upon the absence or presence of vorticity 7 . 

The vorticity within the jet flow will consist of streamwise components 
as well as peripheral components, in general. The surface bounding 
the two regions Is not a streamtube because the turbulent mixing 

action of the jet flow results In entrainment of air from the aero- 
dynamic field. 

Consider the simply connected region t containing a jet. The 
velocity field U is uniquely defined by its divergence and curl 
throughout the region and the normal component of U on the surface 
bounding the region. For most purposes In aerodynamics, the field 
may be presumed to be solenoldal (incompressible)*. This solenoidal 

-y 

field U, as previously shown (Equation (15)), may then bo expressed 

T ; 

^^ e i^ n ^ ire I s h ere assumed to be Incompressible, which 

would be a reasonable approximation for low-speed jets. For near 
sonic or supersonic jets, this would be clearly untenable. 
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as the sum of the two component fields 


u - U S, + E P 


( 15 ) 


The field U St depends on the distribution of vorticity throughout the 

region. However, the vorticity is confined to the jet flow. Thus, 

U^ ( is solely dependent upon the sources within the jet flow. The 

remaining component field depends only on the distribution of the 

-> 

normal component of U on the surface of the region. 

It is possible to determine the aerodynamic field in two ways. 

The first makes use of only the scalar potential, and the second 
utilizes both a scalar and a vector potential. In the following 
discussion, the symbol <f> will be employed to denote the scalar potential 
of the aerodynamic flow field when the field is completely determined 
by only the scalar potential. The symbol i|> will denote the scalar 
potential for a flow field characterized by both a scalar and a vector 
potential. In the latter situation, the scalar potential ij> alone does 
not entirely characterize the flow field. 


A. Method based on the S calar Potential 
Consider a domain chosen in such a way as to exclude the jet flow. 
Refer to Fig. 6. The boundary (I) of the domain then contains the 
surface of the jet flow as one of its parts. Since the aerodynamic 
field within the domain is entirely irrotational and solenoidal, it 
may be completely represented by the potential field component E of 

r 

Equation (10). Recalling Equation (14), the relation for E p is 

Ep = - I VG (X. X,) U, • dS, (14) 

•'s 


I 
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FIGURE 6. DOMAINS FOR AERODYNAMIC FIELD 
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Any solution for II in the subject domain requires specification 

-v 

of the normal component . dS^ on the surface of the domain. For 

most applications, such values may be obtained for the parts of the 

surface not bounding the jet flow. For those parts of the surface 

-> -+ 

bounding the jet, values for • dS-j are available only with a rather 
high degree of uncertainty, even for simple jet configurations. That 

O 

is true of data obtained either analytically or experimentally . 
Furthermore, even if the boundary conditions on the surface of the jet 
could be obtained with acceptable accuracy, the solution of the potential 
problem (either by evaluating Green's function for the domain, or by 
alternate methods) remains complex except for problems with simple 
symmetry. 


B. M etho d Ba sed on t he Stokes Pot entials 
Consider a boundary (II) taken around the region containing 
both the jet flow and the aerodynamic field (Fig. 6). This defines 
a domain for which the velocity field U is represented by + E 
since the vorticity within the jet flow is not zero. This choice 

for the domain eliminates the requirement for knowing the normal 

► 

component of II along the surface of the jet flow. In other words, 
the surface of the jet flow is no longer a boundary of the domain. 
The field in the region exterior to the jet flow is i rrotational . 


Thus, the volume integral for the determination of U r ^ in Equation (15) 
is taken only over the region occupied by the jet flow. 

The formulation of this problem is somewhat more involved than 
the method of the scalar potential. Rewriting Equation (15), the 
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relationship for IJ is 

5 = i J < v i * V x j|pr d, i - hf ’ G <*• V • <*1 <=*> 

■> 

Mow, the solution for involves the vorticity distribution 
within the jet region. Certainly the distribution is not known with 
precision for any turbulent jet. Nonetheless, the integral is the 
value required rather than the vorticity distribution proper, so that 
the uncertainty of the integrand may not be a serious problem in some 
circumstances. That is, if the integrand is only approximated, it 
may still be possible to obtain accurate values for the integral. 

This will be considered in Section IV. As far as the evaluation of 
the Green's function for the domain is concerned, there is more flexi- 
bility compared with the method in the preceding subsection. The 
choice of the domain can be made in any convenient way without requir- 
ing the boundary of the jet flow to be a part of the boundary of the 
domain. 


C_. Vorticity in the Jet Flow 

The vorticity distribution in a turbulent jet is strictly a time- 
varying function. As noted in Sec. II, however, when only the time 
averaged vector field is to he determined, specification of the time 
averaged vorticity field and time averaged boundary conditions are 
sufficient to determine the solenoidal field (Equation (15)). This 
is due to the linearity of the vector field equations. 

For most jets the principal time-averaged vorticity component 
is peripheral, or normal to the streamwise direction of the jet flow. 
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That is, the peripheral component of vortex lines form closed curves 
about the axis of the jet. For a jet in a cross-flow, there will be 
in addition, to this peripheral component, components of vorticity 
in the streamwise direction along the jet. The vorticity lines, re- 
sulting from the peripheral and streamwise components, will then be of a 
helical character. It is advantageous, for the purposes of discussion 
and for computations as well, to consider the vorticity field as 
separated into the sum of peripheral components and a streamwise 
component. Thus, for the vorticity field , 

ft, = u + n (25) 

1 p z v ' 

*f 

where sip is the peripheral component and lies in a surface normal to 
the streamwise direction, and is the streamwise component. 

The magnitude of the peripheral component of vorticity in the jet 
flow has a distribution such as that illustrated in Fig. 7. The 
vorticity on the axis of the jet is zero, and the vorticity approaches 
zero asymptotically near the edge of the jet. Maximum values of 
vorticity are reached at positions between the axis and the edge of 
the jet flow. 

Computations of the aerodynamic field induced by a jet by means 
of Equation (15) requires specification of the vorticity throughout 
the jet flow region. As mentioned before, precise data for the vorti- 
city are not available, either from experimental data, or from analyt- 
ical relations, even for jet flows of the simplest symmetry. However, 
there is not actually a need for precise vorticity data to obtain 
reasonably accurate data for the aerodynamic field, because only an 
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integral of the vorticity is involved iri the computations. It is then 
possible to obtain acceptable approximations for the computation of 
the aerodynamic field if only certain integral properties of the vorti- 
city distribution are known, that is, providing the field points of 
interest are not close to the surface of the jet. At field points 
extremely close to the jet, of course, more accurate data become 
necessary. 

When field points are sufficiently far from the jet flow region, 

* 

the value of | R| in Equation (12) does not vary much over a cross 
section of the jet flow. The integral, then, can be reduced in an 
approximate way to a simpler relation involving only the integral 
of the vorticity and itsrrioments over finite volume elements of the 
jet flow. This approximation leads to useful computational simplifi- 
cations. For some problems, the integral of the vorticity for the 
region can be related to data for the jet flow that is known with 
considerable accuracy. Aerodynamic computations for such jet flows 
employing these methods might be expected to be of comparable or 
higher accuracy than computations previously made with other methods. 

To illustrate the foregoing remarks, consider the contribution 
to the velocity at a field point X due to the vorticity within a small 
section of the jet flow. Refer to Fig. 8. The dimensions of the 
section of the jet flow under consideration are presumed to be small 
in comparison with the distance from the element to the field point. 

The section is a plane slice of the jet flow field of infinitesimal 
thickness. The extension to finite elements of the jet flow will be 
considered later. 
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Equation (12) may be expressed In the form 


“l x 7 1 dT l 


u s, if 

Expanding 1/R as a Taylor's series in X ] , yields 6 ^ 

= X ' x li <r> + ? x 11 x lj iVf 5 Xj ( T> ' ' ' 

~¥~ 

where X = |X| 

The subscripts denote components and the summation convention is 
implied. If only the first and second terms are retained, then in 
vector form 


(26) 


1 

-v 

IX-X, 


1 v J 

J ' X 1 * 7 Y 


Substituting for 1/R in Equation (26), yields the approximate equation 


f a, x 

! - ! j 

7 x - X 1 * v < 7 P 

K 



dx- 


Since X is not a function of the source coordinates, X 1 , this equation 
can be expressed as 


u st “ h 7 x x J Q i dT i + h jf “1 X [ x i • v < 7 ft 

However, the integral of the first term can be expressed as 

Jf dr, = J i>„ dx, - • (X,, «,) dx, 


dxi 


(27) 
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From the divergence theorem, then 



( ^ i -j ^ i ) dx i = 




(28) 


Considering Equation (28) and the expression for the vorticity 
in Equation (25), yields for the integrand of Equation (28) 





+ ^ 2 ) • d$j 


(29) 


It can be concluded that since a is everywhere tangent to the surface 
S of the element then the contribution from the peripheral component 
of vorticity in Equation (29) is zero. If the streamwise component 
of vorticity, « z , is non-zero, then the first term of Equation (27) 
yields 


1 

47 


tt V X X j iQ l dT l " 4ir V X X Jf 


“z dT l 


(30) 


Since the first term yields no contribution of the peripheral 
> 

component to U<^ , the second term of Equation (27) must be used to 
evaluate its contribution. Likewise, if the integral of the stream- 
wise component in the first term is also zero, then the higher-order 
approximation must be obtained, using the second term of Equation (27). 

Consider now the case where the vorticity is entirely peripheral, 
a z =. 0. The first integral in Equation (27) is zero. The expression 
for U<^ in Equation (27), using the vector identity 


v x (<f)A) = <f>v x A + (v<j>) x A 



3 ? 


then becomes 





n l X 1 


7 X dT l 


(31) 


Rewriting the dyadic In the Integrand in terms of skew-symmetric 
and symmetric parts, yields for the Integral In Equation (31) 

J V i ' 7 r d, i * \ f ( Vi • 7 1 - x, • 5 [) <it, 

T ~T 

(“l X 1 7 5T + x i a i * 7 X") dT i 



Noting that the second integral on the right can be rewritten as 


7 ‘ JT («1 *, 


+ X 1 V d, i > • ’r 


and since 


J ^ (X li 8 lj + X lj n ii } dT l = jf 


V * - X 11 X lj V dT l = 0 


only the first Integral (skew- symmetric part) contributes to U r In 

ox 

this case. Using the vector identity 


B (C-A) - A (C-B) = - C x (A x B) 
the skew-symmetric integral can be written as 


/ 


\ ( v J-) x I (X 1 x Qj ) dt-| 


Substituting this value for the Integral In Equation (31) and using 
the vector identities 
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Vx (A x B) = B • VA - A • Vli + AV-B - BV-A 


-> -+ 


v ( A* B ) = A-vB + B-vA + Ax (vx B) + Bx (v x A) 


finally yields , 

* 1 


u, 


St Eir" 


<'T> 


f- - 1 

I x fi-j dx i 


(32) 




com- 


Equation (32) shows that for the conditions considered, that of peri- 
pheral vorticity only, the significant term remaining involves only 
the integral of the moment x u. , over the volume element of the jet. 
It is possible to relate the latter integral to the circulation about 
a vortex ring or an element of a vortex sheet, and then to the velocity 
in the jet flow. If one considers the more general case for v 
f (*i x u.|) d^ jit is possible to relate this term to either a 
bi nation of a source doublet and a vortex doublet. 

For computational purposes, it is possible to replace the volume 
distribution of vorticity in the element of the jet by a discrete, 
closed line vortex of appropriate properties. The properties of the 
line vortex corresponding to the volume distribution of vorticity may 
be evaluated as follows. 

For a plane closed curve C, the line integral 


% =1 
H 2 




X x dX 
1 


1 


(33) 


yields a vector normal to the surface of a magnitude equal to the area 

enclosed by the curve C (Fig. 9). The circulation, r, is defined as 

-> 

the line integral of U 1 around a closed curve C bounding an arbitrary 
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surface in the flow. Thus 


f 


dLi 


• r 

where dL^ is an element of the curve. By employing Stoke's theorem the 
relation between circulation and vorticity may be expressed as 


/ 


dS-j = 


f 


dL 1 = r 


(34) 


where the surface integral is over the surface hounded by the curve C. 

Considering one-half times the integral of Equation (32) and using 
Equations (33) and (34), yields for a closed ring vortex in the plane 
normal to the axis of the jet. 



(X 1 x dX -j ) n p dS 1 


where the surface element dS 1 is that of the cross section of the vortex 

-v 

ring and dX 1 denotes the differential line element along the vortex ring. 
The above expression may be written 


1 


/ 


(X ] x dX-j ) I u 


l 


n dS-| = r A-j 


where « p is the vorticity of the closed ring vortex In the direction 
dX-j and A 1 is the area enclosed by the ring. By considering a series 
of vortex rings along the axis of the jet, an element of a vortex sheet 
is obtained. 


Considering finally v 


r 


v 1 


/ 


► -> 

(X i x q-j ) dfj 


and letting the 


integral term be represented by m, yields 
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v ( V Y‘ * 'i') = (ni * v) ^3 = X [m • (v --y )] + -y (m • v) X 

X XX 

and expanding the above, results in an expression for the field of a 
dipole singulari ty. 



Thus, Equation ( 32 ) can also be expressed as a velocity field induced 
by a dipole of strength m equal to the integral of the moment X-j x iz, . 

The streamwise component of vorticity o is more difficult to 
treat analytically, but the following might be noted. Considering the 
first order approximation expressed by Equation (30), this yields for 
an equivalent line vortex 


1 

47 


v 





1 


where dX^ is a vector element of the line vortex. When r is zero (no 

■ ► 

net rotation of the jet), the second term involving u z in Equation (27) 
must be considered. Liy analogy with the foregoing for it may be 
anticipated that the contribution will be equivalent to that of a 
vortex doublet. 

All of the above approximations are for large R. As previously 
stated, more accurate approximations may be required for points of 
the aerodynamic field close to the jet flow. The foregoing approxi- 
mations provide insight into the nature of the analysis, as will be 
shown in the following sections. 
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IV. THE AX I SYMMETRIC JET IN A QUIESCENT ATMOSPHERE 

As an application of the methods discussed In the foregoing sections, 
the most elementary situation with regard to jet flows is that of the 
axisymmetrlc jet exhausting into a quiescent environment. The vorticity 
in a jot of that type is entirely peripheral. No streamwlse vorticity 
is present, except in the turbulent motion of the fluid. Since only 
the time averaged properties of the flow field are of significance In 
computations of the steady aerodynamic field, turbulent fluctuations 
do not enter the computations. Before examining the computational 
aspects of determining the induced aerodynamic field, the nature of 
an axlsymmetric jet is briefly reviewed below. A synopsis of some of 
the literature pertinent to such jets is included in Appendix A. 

A turbulent jet exhausting into a quiescent atmosphere is generally 
considered to comprise three regions (Fig. 10). Region I, the potential 
core region. Is characterized by a central region of uniform "potential 
flow,"* separated from the aerodynamic field by a turbulent mixing layer. 
The mixing layer grows in thickness downstream from the lip of the 
nozzle. At the end of the potential core region the Inner boundary of 


Such a characterization presumes the dominant flow from the nozzle Is 
essentially uniform. There will always be boundary layers present 
adjacent to the wall. In some situations, the flow from the nozzle 
might have the appearance of a fully developed turbulent pipe flow. 
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FIGURE 10. AXISYMMETRIC TURBULENT JET 
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the mixing layer reaches the axis. The growth of the mixing layer is 
observed to be approximately linear and the velocity profiles, as a 
function of a radial difference normalized with respect to the width 
of the mixing layer, exhibit similarity. At some point downstream 
from the end of Region I, the flow becomes developed in the sense that 
the radial profiles of axial velocity exhibit similarity (Region III). 

That is, the radial distribution of velocity, normalized with respect 
to the velocity at the centerline of the jet becomes a function of 
r/r^ only. The radius r^ is the radius at which the velocity is half 
that at the axis of the jet. In Region III, r^ Increases linearly with 
distance in the downstream direction. The velocity at the centerline 
decreases as 1/Z where Z is the streamwise coordinate along the jet. 

Region II, between Regions I and III, termed the transition region, 
has been found to exhibit non-similar velocity profiles. The decay 
of the centerline velocity begins in the transition region. For 
most purposes, the transition region may be considered to be rather 
short, on the order of 2 nozzle diameters in length. 

The boundary of the jet flow, Introduced in Sec. Ill, was defined 
to enclose the region containing significant time-averaged vorticity. 

The instantaneous position of a jet surface has the appearance shown 
in Fig. 10. The jet surface separates regions of irrotational fluid 
from those of turbulent, rotational fluid. The above definition of 
the jet boundary, then, would be a surface enclosing the extreme positions 
of the Instantaneous surface of the jet flow. 

The mass rate of flow, associated with the jet flow. Increases in 
the downstream direction. Thus, the boundary of the jet flow is not 
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a streamtube since fluid from the environment crosses the boundary to 
become entrained In the jet flow. In the developed flow region of an 
axisymmetric jet, the rate of increase of the mass flow rate per unit 
length along the jet becomes constant. This rate of Increase in the 
mass rate of flow in the jet per unit length of the jet is termed the 
mass entrainment rate, E m< The mass entrainment rate is simply related 
to the normal component of velocity at points on the boundary of the 
jet. For an axisymmetric jet, the mass entrainment rate may be expressed 
as 

E m = 2 " r o U n P 

where r Q is the radius of the boundary and 1)^ is the component of 
velocity normal to the boundary and directed Inward toward the jet 
axis. Recalling that the boundary conditions required for determining 
the aerodynamic field by means of the scalar potential (Sec. 1 1 1 A) in- 
volved specification of the normal component of velocity on the bound- 
ary of the jet flow. It is apparent that those boundary conditions 
are equivalent to specification of the E m . 

Now, as noted in Sec. 1 1 1 A , E m is known for jet flows only with 
a rather high degree of uncertainty. That Is true from both analytical 
and experimental viewpoints. For example, Ricou and Spalding^ examined 
the available data and analytical relations finding that if the entrain- 
ment were expressed as 

E m * K 1 /Rp (35) 

where M is the momentum flux and P the mass density of the jet fluid, 
then values for were found to range from 0.22 to 0.404. Experiments 
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by Ricou and Spalding^ to obtain the value of the entrainment rate 
yielded a value of 0.282 for the mixing of an air jet with air. On 
the other hand Wygnanski^ has found the higher value for corres- 
ponding to the Schllchtlng profile, of 0.404 yielded agreement with 
data he obtained. Other studies appear to show similar uncertainties 
in the entrainment rate. Such uncertainties might be expected because 
of the sensitivity of the entrainment to regions of the jet near its 
edge, where measurements are difficult to obtain with accuracy and 
where theoretical relations there may not be expected to be accurate, 
that is, at least in the boundary layer approximation. 

From the foregoing remarks. It is apparent that however accurately 
the geometry of the problem of computing the aerodynamic field may be 
treated, the boundary conditions are known only with a high degree of 
uncertainty for methods based on the use of the scalar potential alone. 
The following discussion, then, is concerned with similar computations 
based on the use of the Stokes potentials. The use of the Stokes 
potentials for the computation Involves data for the jet flow known 
with higher accuracy. Thus, its use might be expected to yield more 
accurate results under some conditions. 

A. Formulation of the Problem 

Consider an axi symmetric. Incompressible, turbulent jet exhausting 

normally through an infinite plane baffle plate into a quiescent 
■> 

atmosphere. U Is uniquely defined In the semi -infinite, simply connected 

-► 

region t by specifying the curl of U throughout t and the normal compo- 

*-► 

nent of U on the surface S bounding the region t. Since the time- 
averaged flow is assumed to be axlsymmetrlc and without swirling action. 
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the vorticity has only one component u (cylindrical coordinates), as 
shown in Sec. I ID. The vorticity is solely peripheral. This implies 
that the aerodynamic field is also axi symmetric with only r and z 
components. The entire domain of vorticity will be assumed to be that 
of the mixing region around the potential core and of the fully 
developed region of the jet flow. The transition region of the jet 
flow is not considered separately in this formulation. That is, the 
fully developed part of the jet flow is simply presumed to start at 
the end of the potential core. 

As discussed in Sec. II the aerodynamic velocity field may be 
separated into the sum of two fields and E . The boundary condi- 
tions to be satisfied for the domain are that the normal velocity 
component on the plane baffle plate be zero; it will be presumed that 
the normal velocity component at points of the plane within the cir- 
cumference of the jet orifice Is uniform and equal to U., the jet 

J 

velocity. The velocity at distant points approaching Infinity must, 
of course, approach zero. The simplest means for satisfying those 
boundary conditions for this problem is to use the method of images. 

An image of the jet flow is placed along the negative Z axis. Con- 
tributions to the normal component of velocity at the baffle plane due 

-f 

to U $t for the sum of the jet flow and Its image thus vanish. The 

normal boundary condition remains to be satisfied by the E field*. 

P 

*Strictly speaking, to be consistent with the discussion of Sec. Ill, one 
should regard Ep as the sum of the contributions of the Image of the re- 
gion of vorticity plus the Ep field actually considered, for the semi- 
infinite domain. There Is no fundamental difference involved, however, 
so there should be no confusion In the matter. 
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■> 

The normal component of L p is zero at points in the baffle plane outside 
the circumference of the jet orifice, and U. at points In the plane 

J 

within the circumference of the jet orifice. 

The aerodynamic velocity field U for the axisymmetrlc jet has two 
components (r and z) in cylindrical coordinates. Recalling Equations 
(12) and (15) the contributions to the r and z components of U due to 
the vorticity in the jet flow may be expressed as 



U-z,) 

i r i dr i dz i de i 


(r^-r cos o, ) 

r 3 r l dr l d2 i do i 


(36) 

(37) 


where the notation corresponds with that for cylindrical coordinates 
discussed in Sec. III). The nonsubscripted variables refer to the field 
point X, and the subscript 1 denotes values evaluated at the source 
point Xj. The geometry associated with the Equations (36) and (37) is 
Illustrated in Fig. 11. A small segment of the jet is shown in Fig. 11 
as a truncated cone. The integrals of Equations (36) and (37) are 
taken over the points within the jet flow and Its Image. The value R 
in Equations (36) and (37) indicates the magnitude of R = X - X-j , as 
shown in Fig. 11. 


The irrotational and solenoidal part of U, denoted by E p , satisfies 
the boundary conditions for the semi -infinite domain and is represented 
by Equation (14). For the present application, the equations may be 
written in component form as follows. 


E = ^ 
pr 2it 


I. 


(r-r^ cos e-j ) 


r^ dr, de^ 


(38) 
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PZ 


v 


f. 


( 2 - 


Z 1 } 


r, dr, d0 ] 


(39) 


where the notation corresponds to that previously discussed. The 

geometry for Equations (38) and (39) is shown in Fig. 12. The surface S 

surrounding the domain consists of an infinite plane at z = 0 and the 

remaining surrounding surface at infinity. The normal component of 

velocity along the boundary is zero except at the orifice. At points 

in the plane within the orifice circumference the normal component of 

velocity Is presumed to be uniform with the magnitude U., the exit 

J 

velocity of the jet. Thus, the only surface contributing to the surface 
integral is a flat disk of radius r., the exit radius of the jet orifice. 

J 

Integrating over the 0, direction and recalling that fl 0 , Is Indepen- 
dent of 0, , the above equations can be written as follows. 


Sit 
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(z-z,) 
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r,) 2 + (z-z, ) 2 


^r 2 + r 


(r-r, ) 
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+ (z-z,) 2 
+ (z-Z-| ) 2 


E (k) - K (k) 


dr, dz, 


(40) 
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E(k) - K(k) 


dr, dz, (41) 
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"P 2. 
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1M 


7-- 9 dr, 

/(r + r^)2 + z 2 [(r-^) 2 + z 2 ] 1 


(43) 


where K (k) is a complete elliptic Integral of the first kind and E (k) 
is a complete elliptic Integral of the second kind. Derivations of 
Equations (40), (41), (42), and (43) are presented in Appendix B. Note 
that the integrals in Equations ( 40 ) through (43) involve integration 
only over the r^ and z ] coordinates in the jet flow. 


B. Numerical Method 

The integrals in Equations (40) through (43) may be evaluated 
numerically. For the purpose of numerical integration, it is convenient 
to transform the variables of integration in such a way that the domains 
of integration become rectangular. This may be accomplished in the 
following manner. Denote the Inner radius of the mixing layer in the 
potential core region by r^ and the outer radius of the jet by r . 

The geometry for the jet flow is shown In Fig. 13a. 

Introduce the transformation for the coordinates (r^, z 1 ) in the 
mixing layer of the potential core region to the transformed variables 

(x r 


X 1 = (r l" r Il )/(r ol" r Ii ) 
*1 = Z l /r j 

where r n = (z } ). 


Y (44) 


r ol = r o ( z l>* and 
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fj = the radius of the jet orifice 

For the developed region of the jet, the transformation from 
(r-j , z^ to (x-j , ) employed Is as follows. 


X 1 s r l /r ol 
y l = z l /r j 


Y (45) 


The notation corresponds to that described above. The coordinates for 
the surface integration over the orifice area in the plane z 1 = 0, let 


y, - z, = 0 


> (46) 


Finally, for the field point coordinates, introduce the transformation 
from (r, z) to (x, y) as follows. 


x = r/r Q (z) 
y - z/r d 


^ (47) 


The transformations (44), (4b), (46), and (47) transform the original 
domain in the (r, z) plane to the rectangular domain in the (x, y) plane 
as shown in Fig. 13. Region 1 Is the mixing region surrounding the core. 
The fully developed jet flow is represented by region 2. It should be 
noted that the region of integration In the x-j direction extends only 
from 0 to + «.. The region of Integration in the y-j direction contains, 
in addition to the interval (0, + «), the region of Integration over the 


image from 0 to -®. 
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It is convenient. In addition, to normalize the velocity components 
with respect to IK, and the vorticity component with respect to U./r.. 
Thus, the appropriate relations are 


U x ■ U S,/ U j 

E x ■ V U j 

u y 3 u Stz/ u j 

E y * E pz /U j 

'"» r j /u j 

' ■ V u j 


Substituting the relations (44) to (48) into Equations (40) through 
(43) yields the following equations. 


U = 
x 


-^'£ 77 = 

•'-» y 0 / -j r Q x + 


< r 0 r r n > a <y,-y) 


[(r or r Il> X 1 + r n ]} 2 + ^ (y-y,) 

r o * 2 * t(r or r n> x i * r i )] 2 4 ( y - y ,) 2 


| r o x -t< r 0 r r ii * x i 4 r „]} 2 4 r j (y-y,) 

r 0 , “ (y,-y) 


2 ^ (k-j ) ■ K ( k-j ) 


dx 1 dy 1 


r j r - pi 'ol “ vjr 

rrTl 1 / =■ - 

J— J 0/ ( r x * r x ,) 2 + 


o x + r oi x i ) 4 r j (y-y,)' 


I .2 ..2 


r o x ~ 4 r o x i 4 r , (y-y , ) 2 

Z ~2 2 “T E ( kp ) - K (k~) 

(r o x " r o, v 4 r , (y-y,) 2 2 


dx 1 dy ] 


(49)’ 


*imLl n fo 9 tft t10 l 1 S t 5° u e P® r J ormed °ver the jet flow (0 to -) and Its 
ge (0 to -«>). it should be noted that y ^ is replaced by -y, and n 
Is replaced by ~u for the Image Integration] ^ 
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-tffr- 

■'-» *0 /< r x + 
o 

V. 


(r or r i> “ 


t< r or r n> * r n 1 } 2 * r j <y-J'i 


2 2 


r o x ‘ - t(r or r n> x i * r n] 2+ r j (y-yi ) 


{v • [ (r oi' r n ) X 1 + r n f 7 J(y-y }) 


2 E (k,) - K (k,) 


O ’ r o1 0 

o/( r 0 x + r 1 x, ) 2 ♦ rj 


2 (y-y ,) 2 


dx 1 dy 1 


gjj ♦ r j <^i> 2 

(V- r dl x l> 2 + r j <^l> 


? E (k 2 ) - K (k 2 ) 


dx-j dy-j 


E=A f 1 

X Trr n X 1 / o 

Jo /r Q x + r^.x-j r + 


2 2 
r j y 


2 2 2 2 2 2 
V ” r -i x i~ r ’-i y 


’■ITJ-J. 


--jf 1 2x » 

* oyfr x + r ,x. + 


(r o x - r j x l) 2 + r j * 2 


E (k 3 ) + K (k 3 ) 


E (k.) 


(50)’ 


dx 1 (51 ) 


dx-j (52) 


0/(r o x + r^) 2 + y 2 [( r 0 x " r j x l^ 2 + r j 

where for Equations ( 49 ) and (50) the first Integrals on the right side 
are taken over the core region and the second Integrals are taken over 


*The Integration Is to be performed over the jet flow (0 to ») and Its 
Image (0 to -«>). It should be noted that y. Is replaced by - y, and a 
Is replaced by -a for the Image integration. 
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the developed region. The moduli k, , 4 and k 3 for the complete elliptic 
integral functions are 


kf = 


_ 4 _y C (r oi- r n> x i + r n ] 

Il> X 1 + r n Jj 2 + ^ (y-y,) ; 


{v + t< r 0 r r 


(53) 


ko = 


4 r o X r ol X 1 

T'rr 


( V + r oi x i) +r j(y-y]) 


(54) 




4 r o x r ,j x 

X + "j X l^T y 


(55) 


It should be noted that the transformations employed do not strictly 
lead to the integrals in Equations (49) and (50) because the trans- 
formed coordinates are not orthogonal. Nonetheless, the error Involved 
In the approximation of the differential surface element can be shown 
to be of negligible order because of the small spreading angles of the 
jet. Numerical Integration of the equations by this method and by a 
different method not involving a transformation of the domain yielded 
Identical results to four significant figures. 

For the purposes of numerical computation, the numerical algorithms 
for evaluating the complete elliptic integrals were taken from refer- 
ences 12, 13, and 14. Those algorithms were found to be compact and 
of demonstrated accuracy to six significant figures. The algorithms are 
listed In Appendix C. 
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Tht Integrals In Equations (49) and (50) are of the form 



(*!• y^) 4y| 


(56) 


where the variables (x, y) not entering the Integration process are 
suppressed for convenience of notation. The method of Integration 
adopted for the numerical evaluation of the Integrals employed the 
trapezoidal rule for double integrals**. The algorithm may be obtained 
as follows. Let f (xj, y^) be locally approximated by the relation 

f (*!• yj) • « ♦ B yy ♦ t Kj ♦ « *1 y 1 (57) 

In e rectangular element of the domain of dimensions a and b (Refer to 
Mg. 13). Assimm for the purpose of discussion that the origin of 
coordinates Is at the lower left corner of the element and define the 
functions V f V f 2 , and fj as follows. 


f 0 • f (0. 0) f 2 - f (0, b) 

fj • f (a* 0) f 3 ■ f (a, b) 


(58) 


The coefficients In Equation (57) may then be expressed In terms of 


these values. 


« ■ f 0 y * (f r f 0 )/« 

* ■ (f 2 -f 0 ) /b 6 ■ (VVV f o )/ab 


Evaluation of the Integral In Equation (58) then results In the algorithm 


1 • r < f o ♦ f \ ♦ f 2 ♦ V 

where I denotes the numerical approximation for A over the small 
elemental surface considered. 
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The integrals in Equations (51) and (52) involve integration with 
respect to a single variable. The trapezoidal rule was adopted for 
the numerical integration of those Integrals. For an Integral of the 
form 

A ' ’ J 9 ^ dx l 

the numerical algorithm for a small increment of length a may be written 

r * " f (90 + 9,) 

where g Q * g (0) and g-j = g (a), taking the origin arbitrarily at 0 
for the purpose of discussion as before. 

The Integration method employing the trapezoidal rule, either for 
single or double integrals. Is of second order, the error decreasing 
as h where h Is a characteristic size of the element of the domain 15 . 
More precisely, 

A - I = e = c-j a 2 + c 2 b 2 
for the double Integrals, and 

A' - r - £* = Cq a 2 , , 

for the single integrals. This second order behavior of the method has 
been found to hold in the computations made in this study, as evidenced 

by the changes In the solution occurring when the grid size was succes- 
sively halved. 

The division of the different regions Into a network of elemental 
areas was as follows. For the mixing layer of the potential core region 
(first integrals In Equations (49) and (50), a square grid was chosen 
In the Xj , y.j plane. For the fully developed region (second Integrals 
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in Equations (49) and (50), a rectangular grid with aspect ratio of 10 
was chosen, with the longer dimension In the streamwise direction. 
Computations for the fully developed region employing aspect ratios 
of 5 and 2.5 were also made and compared with the aspect ratio of 10. 

These comparisons showed identical results to three significant figures. 

C. Computations for Selected Vortlclty Distributions 
Several theoretical jet velocity profiles were selected for the 
purpose of comparing the aerodynamic fields computed according to the 
foregoing method. For the mixing layer around the potential core, 
the relations suggested by Squire and Trouncer^ (cosine distribution) 
and that of Abramovich v/ were employed. For the developed region 
of the jet, the velocity profiles selected were those of (a) Squire and 
Trouncer^, (b) Schl1chting^ a ind (c) the Gaussian distribution^®. 

The vorticity distributions employed were determined by differentiation of 
radial distributions of the axial velocity. Contributions to the 
vorticity due to the radial velocity component can be shown to be 
negligible.* Computations were made to determine the radial and axial 
velocity components and values for the pressure coefficient, Cp, at 
points along the plane baffle plate and at a series of points in the 
aerodynamic field away from the baffle plate for each of the representa- 
tive velocity profiles. 


*For example, the contribution to the vorticity from the radial velocity 
component of the Schllchtlng profile Is less than 1 percent at _< 
1.2. The contribution increases to 4 percent for 2.8 < r./r^ _< 3.2 a and 
decreases to less than 2 percent at ^/r. = 4.0. 2 
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In addition to the above computations, computations were made 
for the Squire and Trouncer profiles (both that for the potential core 
and that for the fully developed region) to evaluate the sensitivity 
of the aerodynamic field to contributions from various regions of the 
jet flow. For example, the contributions of different axial elements 
of the vortlclty in the jet flow were separately evaluated to deter- 
mine the range of integration In the axial direction necessary to 
obtain accurate values for the velocity components at the field points. 
A separate tabulation of the contributions from different regions In 
the jet flow was made to indicate those regions making the largest 
contributions. 

The contributions to U due to the Irrotational and solenoldal 
part of the field, Ep, were the same for all the jet velocity profiles 
evaluated. A uniform velocity distribution over the area of the jet 
nozzle was employed In determining the E p field. To evaluate the 
Influence of nonuniformities in the velocity profile at the nozzle 
exit plane, several nonuniform profiles chosen to yield the same mass 
rate of flow were examined. 

1. Vortlclty Distributions Corresponding to 
the Selected Velocity Profiles 

a. Mixing Layer of the Potential Core Region . The two velocity 
profiles considered for the mixing layer around the potential core were 

U i 

U zl * z ^ - cos it n] [Squire and Trouncer] (59) 

u z1 = Uj [1 - (1 - n 3/2 ) 2 ] 


[Abramovich] 


(60) 
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where n * (r Q - rj)/^ - rj) and r r < r 1 < r . Evaluating the vorticity 
distributions corresponding to Equations (59) and (60), neglecting the 
contribution to the peripheral component of vorticity due to the term 
3 u r l/3 Zp yields 




01 


IT U 

2 K- r i> 


sin it n 


[Squire and Trounce.r] 


(61) 


and 


fl 0l ■ 


3 U, (1 - n 3/2 ) n 1/2 

"til 


[Abramovich] 


(62) 


For the geometry of the core region, a core length of 12 r^ was 
used as suggested In Ref. 19. For this core length and assuming rj is 
a straight line, an angle of spread of 4.8° for rj was obtained. Apply- 
ing the relation for the conservation of linear momentum to the jet 

between the exit of the nozzle and the end of the core region yields 

>r_ 


n2 „2 _ 9 

u j * r j - 2 • 


/ 


u z, r l *1 


Substituting the value of from Equation (59) Into the above rela- 
tion results In an angle of spread for r Q of 6.7°. 

The results of computations for the radial and axial velocity 
components at a series of points In the aerodynamic field are pre- 
sented In Tables 1 and 2. The values shown in Tables 1 and 2 were 
determined for the vorticity distributions of Equations (61) and (62) 

respectively. The values in Tables 1 and 2 for and 

C C 

are for convenience denoted by U and U respectively. They represent 

* y 

only the parts of the total velocity components (normalized by U^) of 
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TABLE 1. 


y * o 

r o /r j * > 

X 

if 

u c 


X 

y 

l.i 

-7532.0 

0 

2 

-1477.0 

0 

4 

-356.6 

0 

12 

-30.11 

0 

20 

-6.427 

0 

~ — r— — I 

y = 24 

V r j ' 4 - 62 

X 


u c 

i 


FIELD POINT VELOCITIES - CONTRIBUTIONS FROM 

THE MIXING LAYER OF THE POTENTIAL CORE 

[SQUIRE AND TROUNCER (COSINE) VELOCITY DISTRIBUTION] 


y » 12 


6.235 

6.359 

0.2415 

-0.02591 

- 0.01222 


10.57 

4.780 

-0.2443 

-0.09729 

-0.01218 



y = 36 


658.8 

146.3 

20.27 

-0.5179 

-0.1675 


-96.08 

-47.46 

-19.70 

-1.152 

-0.1475 


r o /r j = 6 ' 8 


0.8499 

1.007 

0.4950 

-0.00455 

-0.002530 


1.724 

.9480 

0.01292 

-0.02083 

-0.002620 


y = Z/r. x=r/r o uj = (U^U.) u‘ = (U^/U. 

All tabular velocity ratios have been multiplied by a factor of 10 4 . 







59 


TABLE 2. FIELD POINT VELOCITIES - CONTRIBUTIONS FROM 
THE MIXING LAYER OF THE POTENTIAL CORE REGION 
(ABRAMOVICH VELOCITY DISTRIBUTION)* 


m 

Vj ■ 1 


<< 

It 

ro 

r o /r j a 2 * 42 

X 

U x 

if 


X 

uj 

U C 


X 

y 

" ■ 



X 

y 

1.1 

-7595 

0 


i.i 

797.7 

-127.1 

2 

-1514 

0 


2 

168.1 

-58.01 

4 

-378.1 

0 


4 

22.64 

-22.74 

12 

-33.55 

0 


12 

.5898 

-1.304 

20 

- 7.242 

0 


20 

-.1900 

-.1671 

. 



* » — — • ’ i 

*< 

it 

ro 

-P* 

r 0 /rj = 4.62 


y = 36 

r 0 / r j = 6 - 83 

X 

u c 

u c 


X 




X 

y 



X 

y 

l.i 

7.128 

12.05 


1.1 

.9659 

1.960 

2 

7.256 

5.432 


2 

1.147 

1.073 

4 

2.742 

-.28539 


4 

.5591 

.01210 

12 

-.02946 

-.1103 


12 

-.00521 

-.02344 

20 

-.01386 

-.01382 


20 

-.00285 

-.00294 

*y - z ' r j 

X = 

r/r 0 l 

X o 

II 

u sV u j 

<<>> 

II 

U Stz' U j 


4 

All tabular velocity ratios have been multiplied by a factor of 10 . 
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the aerodynamic field due to the integration over the vortlcity distri- 
bution In the mixing layer of the potential core region. In other 
words, they represent the first of the two Integrals In Equations (49) 
and (50). it can be observed that the magnitude of the radial velocity 
ratio decreases in the radial direction. Similarly, at field positions 
z = 12 r. and greater the magnitude of the axial velocity ratio de- 
creases In the radial direction. The change in the direction of the 
velocity in some cases for r/r Q >_ 4 may be noted. Comparing the com- 
puted results for the vortlcity distribution of Equation (62) [Abramovich] 
with those of Equation (61) [Squire and Trouncer], It can be observed 
that the magnitudes of u£ and were generally higher, some field 
points 11 percent higher, for the vorticity distribution based upon 
Equation (62) [Abramovich]. 

— — The Developed Region. The three velocity profiles considered 
for the developed region were 


= (1 + COS „ ^) 

0 

[Squire and Trouncer] 

(63) 

■ « c EXP [- rf/(2 c| 2^)] 

[Gaussian] 

(64) 


where C 2 = .081 

u 2] * u c /(1 + °* 125 S 2 ) 2 [Schllchtlng] (65) 

where = 18.5 r 1 /z 1 

In the above equations, U c Is the centerline velocity of the jet. Based 
upon a virtual origin of 
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e = 1.2 r. 

J 

where e Is measured from the face of the nozzle exit along the negative 
z 1 axis, the centerline velocity may be expressed by 19 

u c * < u j 13 - 2 + e ) (6*) 

Evaluating the vorticlty distributions corresponding to these 
velocity distributions, again neglecting the contribution to the peri- 
pheral component of vorticlty due to the term 3 U rl /) z 1 ; yields 

U c n r l 

ft ol = " FT s * n v f [Squire and Trouncer] (67) 


a 


el 


r l u 


= — 2 — 1 exp r l^ 2 C 2 z l^ [Gaussian] 


C Z 
c 




171 U c r, 


(1 + 0.125 6 2 ) 3 


[Schllchting] 


( 68 ) 


(69) 


The velocity distributions corresponding to Equations (63) to (65) 
and the vorticlty distributions corresponding to Equations (67) to (69) 
are shown In Fig. 14. The geometry for the region of vorticlty (lateral 
spreading of the jet) for the Squire and Trouncer distribution was 
based upon the selected virtual origin (e * 1.2 r.) and conservation 

J 

of linear momentum for the jet. The angle of spread determined for 
r Q In the developed region (boundary of the jet flow) was 10.4°. 

However, It should be noted that the vorticity distributions for the 
Schlichtlng and Gaussian profiles approach zero asymptotically In the 
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r l /r l/2 


a. VELOCITY DISTRIBUTION 



r i /r |/2 


b. VORTICITY DISTRIBUTION 


FIGURE 14. VELOCITY AND VORTICITY DISTRIBUTION 
IN THE DEVELOPED REGION 
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radial direction. Thus, a large angle of spread, 20°, was selected 
for the jet boundaries of both the Schlichtlng and Gaussian profiles. 

For this angle of spread, the value of vortlclty on the jet boundary 
for the Schlichtlng profile was less than 5% of the maximum vortlclty 
and for the Gaussian profile it was less than 0.1% of the maximum 
vortlcity. Although the jet boundaries differed, the locus of the 
radii (rjJ of the half velocity points was essentially the same for 
all three distributions over the range of Integration. 

The results of computations, using the Squire and Trouncer, Gaussian, 
and Schllchting distributions are presented In Tables 3, 4, and 5 respec- 
tively. The values and represent only those components of the 
total velocity ratio due to the vortlclty distribution In the developed 
region. The values tabulated were determined for the vorticlty distri- 
butions of Equations (67), (68) and (69). The values represent the 
second of the two Integrals In each of Equations (49) and (50). 

It should be noted that the value of r Q employed to normalize the 
radial field position r was based on the lateral spreading of the 
Squire and Trouncer vorticlty region. For this selection of r , the 
radial field positions r = 1.1 and 2.0 fell within and on the nominal 
jet boundaries for the cases of the Schlichtlng or Gaussian profiles In 
the developed region. Thus, calculations for these field positions 
were not made for the Schlichtlng and Gaussian profiles except along 
the baffle plate. The results In Tables 3 to 5 show that the smallest 
contributions to field point velocities are those associated with the 
Squire and Trouncer distribution. The largest contributions come 
from the developed region based on the Schlichtlng profile. It may 




TABLE 3. FIELD POINT VELOCITIES - CONTRIBUTIONS FROM 
THE DEVELOPED REGION [SQUIRE AND TROUNCER 
(COSINE) VELOCITY DISTRIBUTION]* 



All tabular velocity ratios have been multiplied by a factor of 10^. 






TABLE 4. FIELD POINT VELOCITIES - CONTRIBUTIONS FROM 
THE DEVELOPED REGION (GAUSSIAN VELOCITY 
DISTRIBUTIONS)* 


rj r. = 0 


v * 12 


-8.672 

-15.38 

-28.08 

-42.65 


-99.19 

-27.30 


•16.11 


-34.64 


-7.487 -0.6502 


-4.126 


y ■ U ^ 


x = r/r. 


U D = U? 
x U S 




u ; ■ u sx 2 / u j 


All tabular velocity ratios have been multiplied by a factor of 10 4 . 
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‘ r/r o u x ■ u S,r/ u j u 5 * u S, z / u j 

4 

All tabular velocity ratios have been multiplied by a factor of 10 . 
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be also observed that both and uj are negative for all profiles at 
the axial position y * 36. 

2. Regional Contributions to the Aerodynamic Field 

a. Relative Contributions of Different Regions of Vorticlty along 
the Length of the Jet . In contrast with Equations (49) and (50) which 
have limits in the axial Integrations of 0 to «■>, the numerical computa- 
tions for Equations (49) and (50) were limited to a finite region of 
vortlcity. To show the effect of not Including more distant regions 

In the Integrations, the velocity ratios llj| along the baffle plate 
(y = 0) were computed for a sequence of finite regions of vorticlty, 
successively extending the region of Integration farther downstream. 

Table 6 presents the resulting velocity ratios as the region of in- 
tegration was extended In the axial direction (y^ and In the Image 
direction (-y-|). The results shown in Table 6 Indicate that for 
points on the order of 10 diameters aw^y from the jet In the plane 
of the baffle plate, the Integration should extend In the neighbor- 
hood of 100 jet diameters downstream for acceptable numerical accuracy. 

The requirement to extend the Integration In the downstream direction 
diminishes for points closer to the jet, and presumably Increases for 
more distant points. 

b. Contributions from Radial Elements of the Distribution of 
Vorticlty . An Investigation was conducted to determine which radial 
elements (&xj) of the vorticlty distribution made the largest contri- 
bution to the velocity U x at different field positions along the 
baffle plate. The results of this Investigation are shown In Table 7. 
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TABLE 6. FIELD POINT VELOCITIES - EFFECT OF EXTENDING 
THE REGION OF INTEGRATION [SQUIRE AND TROUNCER 
(COSINE) VELOCITY DISTRIBUTION]* 


<< 

II 


Developed Region 

II 

«r» 

X 

u x (y 1 *26) 

U X D (y 1 =65) 

u x (y-|*130) 

u x D (y^eo) 

1.1 

-6.88 

-7.44 

-7.60 

-7.64 

2 

-12.22 

-13.22 

-13.45 

-13.52 

6 

-28.20 

-31.09 

-31.77 

-31.98 

10 

-30.99 

-35.52 

-36.64 

-36.98 

16 

-24.48 

-30.77 

-32.48 

-33.02 

20 

-19.26 

-26.22 

-28.29 

-29.16 


* * m l/r j x * r/r o *1 * Vj u x * U j T / U j U 5 * U S,r /U j 

4 

All tabular velocity ratios have been multiplied by a factor of 10 . 
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TABLE 7. FIELD POINT VELOCITIES - CONTRIBUTIONS 
FROM RADIAL ELEMENTS OF THE JET* 


*< 

If 

o 


Core Region 

r o /r j ■ 1 

AX-j 

aU^(x s 1 . 1 ) 



m 

au£(x=20) 

0-.1 

- 164.6 

- 26.35 

- 4.57 

-0.19 

-0.03 

.1-.2 

- 489.2 

- 81.29 

-15.07 

-0.73 

-0.13 

.2-. 3 

- 781.6 

-135.8 

-27.12 

-1.56 


.3-. 4 

-1010. 

-183.7 

-39.43 

-2.64 

-0.53 

.4-. 5 

-1147. 

-218.3 

-50.19 

-3.81 

-0.79 

.5-. 6 

-1175. 

-233.6 

-57.30 

-4.85 

-1.03 

.6-. 7 

-1084. 

-225.1 

-58.63 

-5.44 

-1.18 

.7-. 8 

- 879.9 

-190.4 

-52.42 

-5.26 

-1.16 

.8-. 9 

- 576.9 

-122.9 

-37.61 

-4.03 

-0.90 

.9-1.0 

- 202.4 

- 47.07 

-14.19 

-1.60 

-0.36 

1 


-1472 

-356.5 

-30.11 

-6.41 


y = 0 


Developed Region 

Vj * 1 

AX-j 

AUJ(X=1.1) 

AU°(x=2) 


AU x (x=12) 

AU x (x=20) 

O-.l 
.1 -.2 
.2-. 3 
.3-. 4 
.4-. 5 
.5-. 6 
.6-. 7 
.7-. 8 
.8-. 9 
.9-1.0 

-0.00472 

-0.05112 

-0.1954 

-0.4628 

-0.8287 

-1.216 

-1.506 

-1.562 

-1.259 

-0.5160 

-00.01 

-00.09 

-00.34 

-00.82 

-01.47 

-02.15 

-02.67 

-02.77 

-02.44 

-00.92 

-00.01 

-00.16 

-00.62 

-01.48 

-02.66 

-03.91 

-04.85 

-05.05 

-04.08 

-01.67 

-00.02 

-00.23 

-00.90 

-02.15 

-03.89 

-05.74 

-07.16 

-07.48 

-06.08 

-02.51 

-00.01 

-00.18 

-00.71 

-01.69 

-03.05 

-04.51 

-05.64 

-05.90 

-04.81 

-01.99 

TOTAL = U° 

= 7.602 

-13.48 

-24.49 

-36.16 

-28.49 


* y = z/r 3 ax, = 4 r,/r 0 x - r/r 0 auj; = aU^/U, aU° - AU° tr /Uj 
All tabular velocity ratios have been multiplied by a factor of 10 4 . 
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0 

The values AIT for the core region and auJ for the developed region re- 
present only that part of the aerodynamic field due to the Individual 
radial elements ax 1 In the respective regions. That Is, AllJj and AU X 
are the results obtained by employing Equations (49) and (50) to In- 
tegrate over the Individual radial elements In the respective regions. 

The total velocity ratios and uj are also presented In Table 7. It 
may be observed that the radial element (ax^ of vortlclty making the 
largest contribution to the field point velocity shifts within a range 
of (.5 - .6) to (.6 - .7) for the core regions. For the developed 
region it may be observed that the radial element making the largest 
contribution Is (aX^) = (.7 - .8). In addition. It may be noted that 
the contribution to the total field point velocity ratios from radial 
element (aX-j) = (.9 - 1.0) varies from 3 to 7 % of the total contribution. 

3. The Potential Field 

The Irrotational and solenoldal part of the velocity field Uj 

denoted by Ep, satisfies the boundary conditions for U as discussed 

in Section III A. The results of computations for the radial and 

axial velocity ratios of E , E and E respectively, are presented 

r * y 

in Table 8. The values E and E shown In the table were determined 

x y 

for a uniform velocity distribution over the area of the jet nozzle 
and represent the results of Integration of Equations (51) and (52). 

It may be noted that the contribution from the potential field to E x 
is significantly lower at y = 12 In comparison with y = 0 and the 
contributions continue to decrease as y Increases. 



TABLE 8. FIELD POINT VELOCITIES - CONTRIBUTIONS 
FROM THE POTENTIAL FIELD* 


y = 0 

r o /r j * 1 


y = 12 

r Q /rj = 2.42 

X 

E x 

E y 


X 

E x 

E y 

i.l 

7396 

0 


1.1 

7.092 

16.09 

2 

1390 

0 


2 

11.08 

13.82 

4 

320.1 

0 


4 

13.17 

8.199 

12 

34.81 

0 


12 

4.691 

.9712 

20 

12.51 

0 


20 

1.957 

.2430 


^1- 

CVJ 

II 

» ! 
l 

Vj ■ 


y = 36 

r 0 / r j = 6.83 

X 

E 

E 


X 

E 

F 


X 

y 



X 

y 

i.l 

1.716 

4.059 


1.1 

.7537 

1.809 

2 

2.711 

3.525 


2 

1.194 

1.577 

4 

3.322 

2.159 


4 

1.480 

.9763 

12 

1.257 

.2723 


12 

.5724 

.1259 

20 

.5312 

.06901 


20 

.2428 

.03204 

* y - z/r 0 

X - 

r/r o 

E x- 

■ V u j 

E y " 1 

: PZ /U J 


All tabular velocity ratios have been multiplied by a factor of 10 4 . 
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To Investigate the effect of a non-uniform velocity distribution 
over the area of the jet nozzle, several non-uniform velocity distri- 
butions were employed. An equation representing velocity distributions 
for turbulent flow In a pipe was selected from reference 17(b) for 
this purpose. 

r i' r l m 

U = U M — -) 
zl m ' r^ ' 


where is the peak centerline velocity. Employing the constraint 
that the mass flow rate be the same as for the uniform velocity profile 


U., yields 

J 


U 4 


U m ■ (nr + 3 m + 2) 


The profiles selected for comparison were based on the following values 


of m and U. 


m’ 


m = .025 
m = .05 
m = .25 


U m = 1.04 U. 
m J 

U m = 1.08 U, 
m j 

U = 1.4 U. 
m j 


Evaluating Equations (51) and (52) with the selected non-uniform pro- 
files, the results obtained were as follows. In comparison with the 
uniform velocity profile, the largest change In field point velocities 
Epr ar, d Ep z resulted from the profile with the largest peak centerline 
velocity of 1.4 Uj representing the largest deviation from the uniform 
profile considered. Even though the field position r/r. =1.1, z/r. = 0 

J J 

showed a 10.4 percent decrease in Ep r> at the field position r/r^ = 2.0, 
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z / r. = 0 the percent decrease in C pr was only 2.6 percent. At axial 
positions {z/r. * 12) the change in field point velocities only 
affected the third significant figure. For the other profiles consider- 
ed, the deviations from uniformity were considerably less, and the 
differences in computed E p values as compared with those for the 
uniform profile were Insignificant. 

4. Summary of Results for the Selected Velocity Profiles 
The combined results for the field point velocity ratios due to 

the mixing layer of the potential core region /U., the developed 

► [) T ^ 

region l'c t /Uj, and the potential flow E p /Uj are presented in Tables 

9 to 11. The field point velocity ratios resulting from a mixing 

layer based upon the relation (Equation (59)), suggested by Squire 

and Trouncer are the same for all combined results. Likewise, a 

uniform velocity distribution over the area of the jet nozzle was 

employed for all results. Thus, ujr/U. and t/U. were the same for all 

combined results. Variations in the combined results may be attributed 

to differences in the vorticity distributions employed for the developed 
region. 

The combined results employing the Squire and Trouncer distribu- 
tion, Equation (66), for the developed region are shown in Table 9, 
where 


The values U x and U y represent those parts of the total velocity ratios 
U r and U, due to the vorticity distribution In both the core and developed 
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TABLE 9. TOTAL FIELD POINT VELOCITIES 
[SQUIRE AND TROUNCER (COSINE) 
VELOCITY DISTRIBUTION]* 


y 

= 0 





r o /r 

" 1 


X 

U v 

E 

U„ 


u. 

E 

u 

c 

u,/u„ 



X 

X 

r 


y 

y 

Z 

p 

t r 

i.i 

-7539. 

7396. 

i 

CO 

♦ 


0 

0 

0 

2.04 

0 

2 

-1487. 

1390 

-97.33 


0 

0 

0 

.948 

0 

4 

-381.1 

320.1 

-61.07 


0 

0 

0 

.373 

0 

12 

-66.52 

34.81 

-31.71 


0 

0 

0 

.101 

0 

20 

-35.34 

12.51 

-22.83 


0 

1 

o 

0 

.0521 

0 


y 

= 12 





r o /r 

. = 2.42 

J 

■ 

U x 

E x 



n 

H 

U z 


S 

U /U 
z r 

i.i 

-180.3 

7.092 

-173.2 


33.98 

16.08 

49.96 


3.25 

-2883 

2 

-114.8 

11.07 

-103.7 


2.600 

13.82 

16.42 


1.10 

-1583 

4 

-64.30 

13.16 

-51.14 


-4.47 

8.199 

3.728 


.263 

-729 

12 

-22.79 

4.690 

-18.10 


-0.8349 

0.9712 

0.1363 


.0328 

-75.3 

20 

-13.09 

1.956 

-11.13 


-0.4279 

0.243 

-0.1845 


.0124 

166 


4 

♦All tabular velocity ratios have been multiplied by a factor of 10 
(Including C p and U 2 /U r ). 
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TABLE 9. (cont.)* 


y 

= 24 



_ . r rAi ■ 4.62 

X 

U x 

E x 

U r 

u y 

E y 

U z 


c p 

U /U 
r r 

1.1 

-127.6 

1.716 

-125.9 

-0.584 

4.059 

3.475 


1.59 

-2760 

2 

-69.95 

2.711 

-67.24 

-1.046 

3.525 

2.479 


.453 

-3686 

4 

-35.07 

3.322 

-31.75 

-2.959 

2.159 

.7997 


.101 

-252 

12 

-11.12 

1.258 

-9.862 

-0.7469 

0.2723 

-0.4746 


.00975 

481 

20 

-5.896 

0.5312 

-5.365 

-0.5583 

0.06901 

-0.4892 


.00290 

912 


y 5 

36 





r o /r 

j = 6.83 

X 

U x 

E x 

u r 

u y 

E y 

U z 


c p 

U /U 
r r 

1.1 

-86.65 

0.7536 

-85.90 

-1.404 

1.809 

.405 


.738 

-47.12 

2 

-47.54 

1.194 

-46.35 

-1.431 

1.576 

.145 


.215 

-31.19 

4 

-23.56 

1.479 

-22.08 

-1.367 

0.9763 

-0.3917 


.0488 

177 

12 

-6.862 

0.572 

-6.289 

-0.9070 

0.1259 

-.7811 


.00402 

1242 

20 

-3.217 

0.243 

-2.974 

1 

-0.5932 

0.03204 

-.5612 


.00092 

1887 


*A11 tabular velocity ratios have been multiplied by a factor of 10 4 
(Including C p and U z /U r ). 
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region of the jet. That Is, they represent the sum of the Integrals 
In Iquatlons (49) and (50). The value C represents the coefficient 
of pressure 


C 


P 


P - P, 
1/2 p U* 


(uj + U*) 


As previously stated, U r and U z are the total velocity ratios due to 
both the vorticlty distribution and the potential field. The combined 
results employing the Gaussian Profile, equation (68), for the devel- 
oped region are shown in Table 10. The values tabulated for Table 11 
correspond to the Schllchting distribution, Equation (69), for the 
developed region. It may be noted that the largest Cp values are 
obtained for the Schllchting velocity profile. In addition. It may 
be observed that the potential field contribution and the vorticlty 
field contribution are of the same magnitude and opposite signs at 
radial field positions 1.1 and 2.0 along the baffle plate. For example, 
examining Table 9, it may be observed that along the baffle plate, 
y = 0, the ratio of E p /U ST varies from .98 at x = 1.1 to .35 at 
x = 20. At y = 12 the ratio varies from a minimum .095 at x = 1.1 to 
a maximum of .24 at x = 4 and then decreases to a value of .12 at 
x = 20. For the field positions at y = 24 and y = 36, the maximum 
ratios .10 and .085 occur respectively at x = 12. 


D. Approximate Computational Models 
In the Interest of making computations of the type discussed In 
the foregoing sections less elaborate for engineering purposes, several 
approximate models of the vorticlty distribution In the jet were 
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TABLE 10. (cont.)* 


y 

= 24 





r 0 /rj =4.62 

X 

U x 

E * 

u 

r 


u y 

E y 


S 

V U r 

4 

-46.62 

3.321 

-43.31 


-0.0734 

2.159 

2.086 

.188 

-482 

12 

-14.01 

1.258 

-12.75 


-0.743 

0.2723 

-0.471 

.0163 

1 

369 

20 

-7.499 

0.5312 

-6.968 


-0.6624 

0.06901 

-0.593 

.00489 

851 


y = 

36 





r o /r i = 6 - 83 

■ 

U x 

E x 

U r 


u y 

E y 

u z 

io 

u /u 

8 


1.479 

-28.65 


-0.9926 

0.9763 

-0.0163 

.0821 

5.93 

9 


0.572 

-8.174 


-1.064 

0.1259 

-0.9381 


1147 

20 


0.243 

-3.885 


-0.740 

0.03204 

1 

* 

o 

CO 

.00151 

1822 


♦All tabular velocity ratios have been multiplied by a factor of 10 4 
(Including C p and U z /U r ) . 
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TABLE 11. TOTAL FIELD POINT VELOCITIES 

(SCHLICHTING VELOCITY DISTRIBUTION)* 


y 

= 0 


r o /r i ■ 1 

X 

u * 

E * 

U r 

u y 

E y 

u 2 

C P 

u /u 

r r 

1.1 

-7543 

7396 

-147 

0 

0 

0 

2.16 

0 

2 

-1497 

1390 

-107 

0 

0 

0 

1.14 

0 

4 

-393.9 

320.1 

-73.8 

0 

0 

0 

.546 

0 

12 

-87.45 

34.81 

-52.64 

0 

0 

0 

.277 

0 

20 

-52.95 

12.51 

-40.44 

0 

0 

0 

.163 

0 


y 

= 12 





r Q /r. - 2.42 

X 

U x 

E 

X 

U r 


“y 

E y 



S 

V U r 

4 

-115.0 

13.16 

-101.8 


7.851 

8.199 

16.05 


10.52 

1576 

12 

-36.98 

4.690 

-32.29 


-0.1705 

0.9712 

0.8007 


0.1043 

248 

20 

-21.59 

1.956 

-19.64 


-0.4626 

0.243 

-0.2196 


.0386 

112 


*AU tabular velocity ratios have been multiplied by a factor of 10 4 
(Including C p and U z /U r ). 
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TABLE 11. (cont.) 


y = 24 r t/ r j * 4*62 


X 

u * 

E x 

U r 

u y 

E y 

U z 


C P 

U /U 
z' r 

4 

-62.02 

3.321 

-58.7 

-0.5849 

2.159 

1.574 


.345 

-268 

12 

-18.59 

1.258 

-17.33 

-0.986 

0.2723 

-0.714 


.030 

412 

20 

-9.962 

0.5312 

-9.431 

-0.869 

0.06901 

-0.800 


.0089 

848 


y 

= 36 





r Q /r. = 6.83 

X 

U 


U 

— 

E 

U 


C 

U /II 


X 


r 


y 

Z 


P 

z r 

4 

-40.0 

1.479 

-38.52 

-1.541 

0.9763 

-0.565 


.1484 

146 

12 

-11.61 

0.572 

-11.04 

-1.396 

0.1259 

-1.274 


.0123 

1150 

20 

-5.501 

0.243 

-5.258 

-0.973 

.03204 

1 

-.941 


.00285 

1789 


4 

*A11 tabular velocity ratios have been multiplied by a factor of 10 

(Including C and U,/U„). 

P 2 r 










81 


examined. For comparison purposes, the models were based on the Squire 
and Trouncer profiles, assumed for the purposes of numerical compari- 
sons to be exact. The radial distributions of vortlclty selected were 
(a) a triangular function, (b) a uniform step-function and (c) a 
delta function. The latter corresponded to replacing the distributed 
vortlclty In the jet by a vortex sheet. 

The motivation for examining the approximate models was to acquire 
some Insight Into the applicability and use of approximate models which 
can reduce the mathematical complexity of Equations (49) and (50). it 
should be noted that the ease of computation Increases as one proceeds 
from the approximate model (a) to the approximate model (c). The cir- 
culation strength per unit length, r', based upon the centerline veloc- 
ity, was chosen to be the same for all approximate models. The tri- 
angular vorticity profile required Integration of a point by point 
distribution of vorticity, which was also required for the Squire and 
Trouncer Model. The second approximate model, the uniform vorticity 
profile, employed an average vortlclty per unit length as opposed to a 
point by point distribution employed In the previous model. In the 
third model examined, the vortlclty was confined to a vortex sheet. 
Thus, the region of Integration, for Equations (49) and (50) employ- 
ing the vortex sheet model was a surface as opposed to a volume. 

1. Triangular Vorticity Profile 

For this model, the sinusoidal distributions of Equations (61) and 
(67) were replaced with the following triangular distributions. 
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V ’ < r r r i )/( V r i> for r iiV r m 


"6, = -"m (V r t )/l V. 1 for r m < r l — r o 


core 

region 


(70) 


= -jl, r,/r for r, < r 
o-j m l m l — m 


u.. = (r -r, )/( r «~0 for r m < r, < r 

0i m o I o m' m 1 — o 


developed 

region 


(71) 


In Equations (70) and (71), r m Is the value of r^ when u 0 = (the 
maximum value of vorticity at an axial position). To determine the 
value of Equation (34) was employed for both the core and developed 
regions. 



(34) 


where S here refers to a meridional section of the jet from the axis 
outward and of unit length along the jet, and C Is a curve bounding 
that section. Assuming that only the centerline velocity contributes 
significantly to the line integral of Equation (34) and substituting 
the values of vorticity from Equations (70) and (71) yields, after 
Integration over r-j 


and 


2 U, 


O « 


Vi 



M <»> 

[ developed"] 
region J 
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2. Uniform Vortlclty Profile 

Employing Equation (34), the average vortlclty In the core region 
and the developed region may be expressed as 



where, as In the case of the triangular model, all contributions to 
the circulation with exception of the centerline velocity are treated 
as negligible. The average vortlclty for the core region Equation (74) 
extends to the edges of the mixing layer. The average vortlclty In 
the developed region Equation (75) extends to the edge of the jet. 

3. Vortex Sheet Model 

Employing Equations (74) and (75), letting (r Q -rj) approach zero 
but maintaining a constant product of n 0 (r Q -rj) yields for the core 
region and developed region 

n 9 1 (V r I ) a " U j 6 ( r m) [core ^glon] (76) 

U e 1 (V r I> = ' u c 6 [developed region] (77) 

The strength of the vortex sheet per unit length In the z 1 direction, 
denoted by y , is then given by Uj or U c> respectively. This conical 
sheet of vortlclty was placed along the half-velocity line of the jet. 
For the core region, the angle of the half-velocity line was taken to 
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be zero. The half-velocity line for the fully developed region was 
chosen to diverge from the virtual origin of the jet at a half angle of 

5°. These choices correspond approximately with those observed In 

20 

turbulent jets. The vorticlty distributions for the core region 
Equations (70), (74), and (76), and for the developed region Equations 
(71), (75), and (77) are shown In Fig. 15. The distributions, 

Equations (61) and (67), obtained from the Squire and Trouncer (cosine) 
velocity distribution are also shown. The ordinate ( r 0 ~ r i)/ U j 
and the abscissa (r Q -r.j )/ ( r 0 “ r l ) are ^ or core re 9'* on ‘ ** should 
be noted that the vorticlty of the vortex sheet (delta function) Is 
Infinite. It may be observed that the triangular profile represents 
the best approximation In comparison with the uniform profile and the 
vortex sheet. 


4. Result of Computations Based on the Approximate Models 
The results of computations, using the triangular profile, the 
uniform profile, the vortex sheet model and the Squire and Trouncer 
(cosine) reference profile are presented In Table 12. The values 
tabulated were determined for the vorticlty distributions of Equations 
(70), (71), (74), (75), (76), (77), (61) and (67). The values U x and 
Uy represent both integrals in Equations (49) and (50) respectively. 

It may be determined from Table 12 that the ratios of the contribution 
of the individual approximate models to that of the Squire and Trouncer 
(cosine) velocity profile have the following trend. For y * 0 the 
ratio uj/u: ef varies from .98 at the point x * 1.1 to .92 at x s 20. 

In contrast, the ratio uf varies from 1.004 at x * 1.1 to 1.55 

A A 



{2»i (r 0 -ri) 
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( r o — r i) /( r 0 — r i) 
r l /r l/2 


FIGURE 15. APPROXIMATE MODELS FOR THE VORTICITY 
\ DISTRIBUTION 
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TABLE 12. FIELD POINT VELOCITIES - CONTRIBUTIONS 
FROM THE COMBINED MIXING AND DEVELOPED 
REGION* 



/ = 0 






Vj * 1 

X 

"J 

U 11 

X 

u m 

X 

u« f . 

X 

4 

u 11 

y 


u; ef. 

1.1 

-7412 

-7573 

-7381 

-7539 

0 

0 

0 

0 

2 

-1456 

-1525 

-1393 

-1487 

0 

0 

0 

0 

4 

-365.1 

-414.9 

-326.3 

-381.1 

0 

0 

0 

0 

12 

-61.43 

-93.53 

-46.71 

-66.52 

0 

0 

0 

0 

20 

-32.67 

-54.81 

-23.33 

-35.34 

0 

0 

0 

0 


y 

= 12 






r 0 /r. = 2.42 

n 


U 11 

X 

u 111 

X 

“I*- 

m 

in 

m 

EH 

1.1 

-160.2 

-520.3 

-94.26 

-180.3 

41.51 

-19.3 

55.63 

33.98 

2 

-102.4 

-265.2 

-62.37 

-114.8 

6.65 

23.33 

14.25 

2.6 

4 

-58.72 

-116.78 

-40.48 

-64.30 

-2.95 

1.99 

-0.36 

-4.47 

12 

-21.26 

-37.10 

-16.06 

-22.79 

-0.6307 

-1.009 

-0.2105 

-0.8349 

20 

-12.29 

-21.06 

-9.454 

-13.09 

-0.3603 

-.8537 

-0.1940 

-0.4279 


*1 = The triangular profile, II * The uniform profile. III * The vortex 
sheet model, ref. = The Squire and Trouncer (cosine) profile. All 
velocity ratios have been multiplied by a factor of 10*. 
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*1 = the triangular profile, II a The uniform profile, 

III * The vortex sheet model, ref. = The Squire and Trouncer (cosine) 
profile. All velocity ratios have been multiplied by a factor of 1CH. 
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UISCUSSION 
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The results of the computations for the axisynmetric jet serve to 
illustrate the character of the method of computation employing the 
Stokes potentials. For example, the relative contribution of the 
potential part of the field E p as compared with the part due to 
vorticity is of some interest. In principle, at least, the potential 
part may be determined with considerable accuracy, for example, for 
more general jet flows, whereas the vorticity contributions are likely 
to be less accurately known. Comparisons of the computed results with 
the available experimental data and with other theories illustrate the 
strengths and weaknesses of the method of computation and of the 
existing theoretical velocity profiles for the jet flow. 


A. Contribution of E Relative to U 

_[> Sjt_ 

The results presented in the previous section. Sec. IV, (Tables 
9 to 11), show that Ep is opposite in sign to U St along the baffle 
plate for the example considered. The field point velocity U along 
the plate is determined by the difference between the magnitudes of 
Ep and U $t . If these values, E p and U st , are of approximately the 
same magnitude, errors in Ep or are magnified in the computed 
resultant U. Errors of this type are most significant near the 
source, and they become insignificant sufficiently far from the 
source. It may be observed from Table 9, for example, that errors of 
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this type are likely to be important at X = 1.1 on the baffle plate. 

At this position, a 2 per cent increase in the magnitude of 
U (U = U /U ) increases the magnitude of U r (U r = radial component 

X X Sty* J 

of U) from 143 to 294, other terms involved being considered constant. 

In contrast to this magnification of error, at field points away from 
the source (Table 9) the relative contribution of E p to the total field 
point velocity rapidly decreases. That is, neglecting the contribution 
from E entirely at y = 36 results in a maximum error of about 1 per 
cent in the field point velocity ratio U/IK. 

B. Comparison of Computed Results with 
Experimental Bat a and Other Theories 

1. Comparison with Experimental Data 
Experimentally determined pressure coefficients on a baffle plate 
in the plane of the nozzle may be compared with the theoretical pre- 
dictions of the coefficient of pressure for the selected vorticity 
distributions. The available experimental data include results 
obtained by Wygnanski 11 and by Gentry and Margason 21 . The corresponding 
experimental and theoretical values of C p are shown in Table 13. A 
comparison of the experimental data, C 2 and C p , with the theoretical 

values for C p are shown in Fig. 16. 

Comparing the experimentally determined pressure coefficients, 
the values determined by Wygnanski, C 1 and C 2 , are lower than those 
obtained by Gentry and Margason, C p and C p , by more than a factor of 
2 at the points close to the source. This disagreement between the 
different experimental data could have been a consequence of differing 
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Table 13. EXPERIMENTAL AND THEORETICAL VALUES OF -C x 10 4 * 

P 


y = o 

EXPERIMENTAL RESULTS 

COMPUTED RESULTS 

V r j = 

1 

x 

c 1 

c 2 

c 3 

c 4 



c 6 

C D 

c 8 



p 

p 

p 

p 


p 

p 

P 

P 


l 



16.3 

9.6 







1.1 






2.04 

2.10 

2.16 

4.43 


2 

2.9 

2.5 

6.1 

4.4 


.948 

1.04 

1.14 

2.09 


4 

1.2 

.8 

2.0 

1.4 


.373 

.417 

.546 

.908 


8 



1.0 

.7 







12 






.101 

.144 

.277 

.302 



*1 Wygnanski's data for Reynold's no. of 53,800 

2 Wygnanski's data for Reynold's no. of 51,500 

3 Gentry's and flargason's data for a plenum- chamber to ambient 
pressure ratio of 1.32 

4 Gentry's and Margason's data for a plenum- chamber to ambient 
pressure ratio of 2.04 

5 Squire and Trouncer velocity profile (Table 9) 

6 Gaussian velocity profile (Table 10) 

7 Schlichting velocity profile (Table 11) 

8 Schlichting velocity profile for developed region and Abramovich 
velocity profile for core region (Tables 5 and 2) 


v 


ORIGINAL PAGE IS 
OF POOR QUALITY 





93 


velocity profiles at the exit orifice. For example, for the data of 

C 

Table 9, a 2 per cent reduction of E would yield a C value of 9.6 

x p 

as opposed to 2.04 at X = 1 . 1 and a value of 1.56 as opposed to .948 
at X = 2.0. From the results of Section IV 3, showing the sensitivity 
of Ep r to non-uniformities in the velocity profile at the nozzle exit, 
it is apparent that such a change could easily be accounted for. 

Changes in theoretical E values become less significant as x increases 

A 

beyond about 2. The different experimental data also tend to show 
better agreement at the more distant field points. It may be observed 
in Table 13 and Fig. 16 that the theoretical magnitudes of Cp are 
generally lower than the experimental values. 

O 

The pressure coefficient C° shows reasonably close agreement with 
Wygnanski's data, as might be anticipated because Wygnanski's data 
had been found to compare well with other results computed using the 
Schlichting profile (Equation (65)) for the developed region of the jet. 
The use of the Abramovich relation (Equation (60)) for the velocity pro- 

g 

file in the mixing layer in the determining of C° should be expected to 
yield more accurate results as compared with the cosine profile of 
Squire and Trouncer (Equation (59)). Experimental data for the velo- 
city profile in the mixing layer have been found to agree quite closely 
with the Abramovich relation . It may be noted that the computed Cp 
magnitudes close to the nozzle are significantly smaller for the 

cosine profile as compared with that of Abramovich, accounting for 

5-7 

part of the difference between the magnitudes of Cp and the 

1-4 

experimental data, Cp 
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On the basis of the above comparison, it may be recommended that 
computations according to the theory making use of the Stokes 
potentials are likely to be reasonably accurate, if 

(1) the velocity profile of the jet at the nozzle exit 
plane is accurately known, 

(2) the Abramovich relation. Equation (60) , is employed for 
the mixing layer of the potential core region, and 

(3) the Schlichting profile, Equation (69) , is employed for 
the fully developed region of the jet, with compu- 
tations involving integration only within a 20° 

half- angle cone. 

These recommendations apply, of course, only to the example con- 
sidered, that of an axi symmetric jet exhausting normally through a 
plane baffle plate. Nonetheless, the sensitivity to those parameters 
observed for this case seem likely to be applicable to more general 
jet flows. 


2. Comparison With Other Theoretical Models 
Wygnanski^, employing a distribution of sinks along the axis of 
symmetry to satisfy approximately the boundary conditions on the 
boundary separating the jet flow from the aerodynamic field, obtained 
the following equation for the pressure coefficient along the baffle 
pi ate . 


9 

P 


— TT 
4irV 


0.715 


0.514 

[1 + ( x /12) 2 ] 1/2 


(78) 


where x = r/r. 

J 
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In developing the model, Wygnanskl chose the value of 0.404 for the 
numerical constant Kj , in the developed region Equation (35), correspond- 
ing to the Schllchting velocity profile. Values for this constant have 

been found to range from 0.220 to 0.404 as noted In Sec. IV. 

22 

Stewart , also examined the Irrotational flow field external to 
a jet. From similarity considerations, Stewart expressed the mean flow 
near the axisymmetric jet boundary as = C/p, where 4>, p are the 
colatitude angle and the radius in a spherical coordinate system 
(Fig. 3) with the virtual origin of the jet located at the origin of 
the coordinates and with the flow along the Z axis. For a jet emerging 
normally through an Infinite plane baffle plate. It was found that the 
velocity potential <p could be expressed by 

<P s C tan <j> 0 In (p sin <J>) (79) 

where <j> Q = the half cone angle containing all of the turbulent 
region of the jet (about 12.5°) 

Employing Equation (79) and Equation (35) with K-j = 0.404, the correspond- 
ing pressure coefficient along the baffle plate (<J> = n/2) may be express- 
ed as 

Cj° = (-0.252/x) 2 (80) 

where x = r/rj 

A comparison of Wygnanskl 's results, the results according to 
Stewart's theory Equation (80), and the results of computations employing 
Stokes potentials with Schllchtlng's velocity profile for the developed 
region Is presented In Table 14. 
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Table 14. COMPARISON WITH OTHER THEORETICAL MODELS - CpXlO 4 * 


<< 

II 

o 




r 0 /r. * 1 

X 


* 

C 9 

P 

4 ° 

1.1 

2.16 

4.43 

8.58 

11.1 

2 

1.14 

2.09 

2.77 

3.34 

4 

0.546 

0.908 

, 

0.755 

0.836 

12 

0.277 

0.302 

0.219 

0.093 

20 

0.163 

0.171 

0.129 

0.033 


*7 


8 

9 

10 


Schlichting velocity profile for developed region 
and Trouncer velocity profile for the core region 
Schlichting velocity profile for developed region 
velocity profile for core region (Tables 5 and 2) 
Wygnanski, Equation (78) 

Stewart, Equation (80) 


and Squire 
(Table 11) 
and Abramovich 
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As mentioned in Sec. VB2, the best agreement with the experimental data 

O n 

may be expected for C as compared with c'. Wygnanski's relation 

9 ^ 

Equation (78) for C^, yields values which had been shown to fit rather well 

11 a 

his experimental data , also corresponds well with C° as might be 

expected. The results obtained from Equation (80) , C 10 , agree fairly well 

g ^ 

with Cp at small values of x but decrease in magnitude more rapidly as 

x increases to larger values. This difference in and is most 

P P 

likely inherent in the different modeling of the jet. Wygnanski 

employed a line of sinks and Stewart's representation corresponded to 

application of the boundary conditions on the surface of the jet. The 

core region was not represented by Stewart's model at all. Rather, 

the developed region was presumed to originate at the point source. 

On the other hand, Wygnanski's analysis does not strictly represent 

the actual geometry, as does the present theory or Stewart's theory. 

8 9 

Therefore, the values of Cp and Cp as might be expected show differences 
at small values of x. 

C. Comparison of Computations for the Approximate Models 
with the Computations for the Reference Model 
(Squire and Trouncer Profile) 

The three approximate models for the radial distribution of 
vorticity in the jet (triangular, uniform, and vortex sheet). Section 
IV D, were chosen to have identical values for the circulation per 
unit length. The common value of r' corresponded to that for the 
Squire and Trouncer velocity profiles for both the potential core and 
developed regions of the jet. The computed values for C p along the 
baffle plate. Table 12, compared quite well with the reference values 
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except for the uniform and vortex sheet models at points distant from 
the axis. As could be anticipated, the results for the triangular 
profile agreed more closely with those for the reference computation 
even at the distant points. The agreement of all of the results at 
points nearer the axis could be expected because 

(a) the C values close to the axis depend strongly on the 

P 

form of the mixing layer in the potential core region, 
and 

(b) the circulation per unit length and the moments of the 
vorticity distributions for the approximate models of 
the potential core layer closely corresponded to that 
for the reference profile 

On the other hand, for computations downstream and at points far 
from the axis, where contributions from the developed region of the 
jet become significant, the agreement between results for the approxi- 
mate models and for the reference profile was significantly poorer, 
except for those of the triangular profile. The disagreement may be 
attributed to the difference between the moments of the vorticity 
distributions for the uniform and vortex sheet models, as compared 
with that of the reference profile. Those differences were greatest 
for the developed region with the geometry chosen for the computations 
With the discussion of Sec. Ill C in mind, the proper choice for the 
radial position of the vortex sheet, for example, would be such that 
the moment of vorticity distribution is maintained in agreement with 
that for the reference profile. The radial position of the vortex 
sheet was simply chosen to coincide with the half-velocity points in 
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the jet which also corresponded with the points of peak vorticity (Fig. 
15). For more accurate results, then, the divergence angle of the cone 
joined by the vortex sheet should have been larger. 

Even though the approximate models yield somewhat different 
results as compared with those for the reference profile, the agreement 
could be acceptable for some computational purposes for all of the 
approximate models. The triangular profile best represents the data 
for the reference profile as could be expected. The vortex sheet model, 
even with the poor choice for its geometry, did yield the reference 
values to a reasonable approximation. Because the computations for 
the vortex sheet model are by far the simplest, the latter comparison 
suggests the use of such a model for most computational purposes. 


D. Evaluation of the Integral Moment of Vorticity 
The vorticity magnitude at a field point has been shown to be 
directly proportional to the magnitude of the integral moment of 
vorticity 


I 


m 



-> 

x dx^ 


(81) 


Recalling Equation (32) , one has 







(32) 


Now, for the examples considered in the previous sections, involving 
an axi symmetric jet, the vorticity is entirely peripheral. I may 
then be determined as a summation of integrals of the form in Equation(81) 
for small slices of the jet as discussed in Section III C. This enables 
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the evaluation of the significance of local contributions of the 
vortlclty along the jet to the field point velocity. 

Considering a small section of the jet flow (refer to Fig. 8) for 
the axi symmetric jet, the magnitudes of an Integral moment of vortlclty 
per unit length along the jet may be expressed as 


I 



s 




(82) 


The small section Is considered to be a plane slice of the axi symmetric 
jet flow field of Infinitesimal thickness. 

Introducing the following transformations Into Equation (82) 


A = a ei < r o ‘ r i )/u j 
“ ■ < r 0 - V (r o ' r i> 


[core region] 


A ■ n el V u j 

“ ■ Vo 


y [developed region] 


yields for the core region 


Vo 


/ 


r 2 (r O " r i ) (r C 

[i 5_ r _L y + _c 


r : )l 2 

] A d y 


(83) 


and for the developed region 


m 


2 " U J r o‘ 


/ 


y A dy 


(84) 
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The values of the integral, Equation (83), v/ere calculated for the 
mixing layer of the potential core region based on the Squire and 
Trouncer profile, Equation (61), and on the approximate models, Equations 
(70), (74), and (76). The values of the integral. Equation (84), were 
calculated for the developed region based on the Squire and Trouncer 
profile Equation (67), the approximate models Equations (71), (75), and 
(77), the Gaussian profile Equation (68), and the Schlichting profile 
Equation (69). The results of these calculations for the mixing region 
were as follows. 


where 


I 

m 




-1 





Y 


y = .297 [Squire and Trouncer profile] 
.290 [Triangular Vortlclty profile] 
.333 [Uniform Vortlcity profile] 
.250 [Vortex Sheet] 


(85) 


For the developed region, the results of these calculations were 


2TTU.ro 


7 = y 


(86) 


where 

Y = -.297 [Squire and Trouncer profile] 
-.290 [Triangular Vortlclty profile] 
-.333 [Uniform Vorticlty profile] 
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-.250 [Vortex Sheet] 

-.390 [Gaussian profile] 

-.500 [Schllchtlng profile] 

It may be observed for the core region that the magnitude of the 
integral moment of vorticlty Equation (85) Is approximately the same for 
all four profiles. Note that the third term on the right side of 
Equation (85), differing slightly for the different profiles, makes the 

I 

smallest contribution to I m . This provides confirmation of the statement 
in the previous subsection C concerning the close correspondence of the 
ys for the different profiles In the core region. 

In contrast to the above results for the core region, the Integral 
moment of vorticlty Equation (86) differs significantly for the six pro- 
files considered. This substantiates the explanation presented in sub- 
section C for the disagreement among the computations for the different 
profiles, where contributions from the developed region of the jet become 
significant. 

Cased upon the proportionality a y as previously discussed, 

the ratio of the Integral moments and the field point velocities for 

the developed region for different profiles should be approximately 

proportional. Selecting the values of U and U for the Gaussian profile 

x y 

Table 10 and the Squire and Trouncer profile Table 9 at y * 36, x = 20 
yields a ratio of 1.28 for ly's. The ratio of ys for the profiles Is 
1.31. Thus, the proportionality Is fairly close. Selecting other pro- 
files similar results are obtained. 

Recall that the computed values of C p as shown In Tables 10 and 11 
were based upon a half-angle spread of 2)° In the developed region 
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of the jet for both the Gaussian and Schlichting profiles. It should 
be noted that a more acceptable half-angle of spread would have been 

o 9 (*>) 

10 to 15 . Most experimental observations, for example, indicate 

the maximum spreading of the jet boundary to lie within that range of 
half-angles*^. 

In order to evaluate the effect of reducing the half-angle to 
correspond more closely with the observations, the half-angle of 10°, 
used for the Squire and Trouncer profile, was also employed for the 
Gaussian and Schlichting profiles. The spreading angle of the half- 
velocity line was left unchanged for the profiles. Reducing the half- 
angle of the spreading of the jet boundary, then, had the effect of 
neglecting the vorticity associated with the Schlichting and Gaussian 
at points outside the boundary. The results of the computations for 

both the Gaussian and Schlichting profiles yielded lower C values 

P 

than those shown in Table 9 for the Squire and Trouncer profile. 

It appears that all three of the profiles considered are not 
truly representative of the observed vorticity distribution in the 
developed region of the jet. The Squire and Trouncer profile with a 
reasonable half-angle of spread of 10° yields values of C p which are 
too low in comparison with experimental data (Table 13). The Gaussian 
profile, even for a spreading half-angle of 20°, yielded results below 
the experimental data. The best agreement with the experimental data 
employing the Schlichting profile was found for a half-angle of spread 
for the jet of 20°. In comparison with typical observed spreading 
angles, a half-angle spread of 20° is not reasonable. However, as 
previously mentioned, even at this large angle of spread the value of 
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vorticity al the boundary of the jet flow for the Schlichting profile 
is still about 5 per cent of the maximum. A more reasonable profile 
with the same general shape for the developed region would be charac- 
terized by a vorticity distribution such that the half-angle of spread 
is about 10 to 15°. In addition, it would be expected that the moment 
of vorticity should be larger, either due to a larger peak vorticity 
or a redistribution of the vorticity profile. 



VI. CONCLUSIONS 


The following conclusions may be drawn regarding the computation 
of the aerodynamic field of a jet according to the method of Stokes 
potentials (scalar and vector potentials). The conclusions, for the 
most part, are supported by the application of the method to the 
axi symmetric jet exhausting through a plane baffle plate. 

(1) The most significant information required for the compu- 
tations are the velocity profile at the exit plane of the 
nozzle, and the integral moment of vorticity along the jet. 

It was found that the potential part of the velocity in the aero- 

-> 

dynamic field E p w as significant for most of the points computed. 

This was especially true for points near the origin of the jet. The 
need for the integral moment of vorticity I was anticipated from 
consideration of the eS/mptotic behavior of the theory in Sec. Ill C. 
This is strictly true for jets for which the only vorticity present is 
peripheral as in the case of the axi symmetric jet. For jets of a more 
complex nature, there would also likely be a requirement for knowledge 
of the integral moment of streamwise vorticity in the jet. For the 
axi symmetric jet, reasonable values for the integral moment of vorticity 
could be obtained from values for the velocity on the axis and the 
radius of the ha If- velocity point as functions of distance along the 
jet, without further knowledge of the jet velocity profile. 
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(?) At field points on the baffle plate within approximately 

one diameter from the edge of the jet, the contribution 

u{j T due to the developed region of the jet was negligible; 

at all other field points considered, to approximately 

20 jet radii from the region of the jet In either the 

streanwise or lateral direction, all three components E , 

-V "d 

U<j T and U St were of comparable magnitudes. At more distant 
field points, the developed region tends to dominate, as 
might be expected. 

Since both E and remained relatively large, but of opposite sign, 
at distances to approximately 6 jet radii from the edge of the jet, 
both must be known with relatively high accuracy for field points on 
the baffle plate within that range. This significance of the velocity 

21 

profile of the jet at the nozzle exit has been observed experimentally . 
Theoretical considerations have been heretofore limited to speculations 
as to the Influence of the jet velocity profile on the mixing layer 
and consequent changes in the entrainment rate. 

(3) The computations tended generally to underestimate the 
Cp magnitudes observed experimentally, except for the 
case where the Schlichting profile was employed, 
together with an overall spreading half-angle of 20° 
for the jet flow, which compared well only with 
Wygnanskl's data. 

The comparison with experimental data indicates a need for improvement, 
both in the experimental data and In the analytical characterization 
of the jet flow. The experimental data available were not consistent, 
differing by more than a factor of 2 (perhaps 4, considering the 
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Ricou and Spalding results). Moreover, the Schlichting profile, 
yielding the most favorable comparison, is not compatible with exper- 
imental observations of the spreading angle of the jet boundary. The 
analytical results were found to be comparable with both the experi- 
mental data and the other theories, within the limit of the state of 
the art. 

Theories relying on the scalar potential yield acceptable and 
simpler relations for the axisymmetric jet. The method of Stokes 
potentials however, offers the added benefit of explicitly accounting 
for nonuniformities in the initial velocity profile of the jet, 
simpler computations for more complex jet flows, the prospect for 
developing simpler vortex sheet models, and less reliance on data 
that are difficult to obtain for the more complex flows (for example, 
reasonable results may be obtained from measured values for U c and 
r, , as opposed to determinations of entrainment required for the 
method employing the scalar potential alone). 

(4) The approximate vortex sheet model was found to yield 
results comparable to those obtained with the complete 
vorticity profiles, and thus may be recommended for 
computational purposes. 

The vortex sheet model might be expected to be useful because (a) it 
should yield precise results asymptotically if I m is properly chosen, 
and (b) for regions of the jet close to the nozzle, I is known with 
some accuracy because the mixing layer is relatively narrow there and 
grows only gradually. The computations for the axisymmetric jet 
supported this view. There appears to be little justification, con- 
sidering the state of the art, to employ more complex models. 
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APPENDIX A 


Annotated Bibliography on Axl symmetric Jets 
Experimental and theoretical Investigations of axlsymmetric 
turbulent jets have been far too numerous to attempt a complete review 
here. Rather, a synopsis of some of the studies Is presented In the 
form of an annotated bibliography. More complete discussions may be 
found in Abramovich 9 , Hinze 19 , Townsend 7 , and in review papers (e.g., 
Halleen 25 , Skifstad 8 ). 

1. Taylor, G. I., "Eddy Motion in the Atmosphere", Phil. Trans. 

Roy. Soc. London , Series A, Vol. 215, 1915, pp. 1 - 26. 


Taylor developed a purely phenomenological theory concerning 

vortex motions. He assumed that the vorticity could be considered 

as a transferable quantity. That is, the vorticity would be conserved 

along the path of a lump of fluid over a certain distance. Based 

upon this assumption, he obtained the relationship 

. 2 


a 


£ 

2 


dU- 


1 


dx. 


where r is the coefficient of eddy diffusion for vorticity, is 

a 

Taylor's mixing length for vorticity, and dU^/dx 2 Is the time mean value 
of vorticity for two-dimensional flow uniform In the x-j direction. Taylor 
extended this theory to three-dimensional flow. This three-dimensional 
theory Is known as the modified vortl city-transport theory. 
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2. Zlmm, W., "Uber die Stronungsvorgange 1m frelen Luft strahl," 
Forschungsarbeiten aus dem Gebiete des Ingenleurwesens . Mo. 234, 

Zlmm conducted an experimental Investigation of an axlsymmetric 
jet issuing into a quiescent atmosphere. In the Zinin experiments 
the jet exit velocity was non-uniform, that is 

U n/ U cp = 1 ‘ 1 


where U m is the peak exit velocity and U cp is the exit velocity close 
to the boundary of the jet. The experimental data for the decay of 
the centerline velocity in the developed region of the jet could be 
approximated by the expression 

U c /Uj = 0.96/ (az/rj) 


where U c is the centerline velocity, Uj is the average jet exit velocity, 

z Is the axial distance downstream from the jet exit, rj is exit radius 

of the jet, and a is a numerical constant with a value of 0.070. 

3. Prandtl , L., "Uber die ausgeblldete Turbulenz," ZAMM, Vol. 5, 

1925, pp. 135-139 (also Proc. Ind. Intern. Cong . Appl . Mech., 

Zurich, 1926, pp. 62-75.] “ 

Prandtl introduced a phenomenological theory concerning momentum. 

He assumed that the momentum could be considered as a transferable 
quantity. That is, the momentum is conserved along the path of a lump 
of fluid over a certain distance. Based upon this assumption he ob- 
tained the two-dimensional relationship 


e m = 


dU 1 

dx^ 
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where c m is the coefficient of eddy diffusion for momentum, is 
Prandtl's mixing length for momentum, and dU^/dx 2 is the difference 
in mean velocity, between the end of the path and the start of the path 
from where the lump originated, divided by the path length. 

4. Tollmein, W., "Berechnung Turlenter Ausbreltungsvorgange," 

ZAMM, Vol. 6, 1926, pp. 468-478 (also NACA TM 1085 [1945]) 

Tollmein published an analytical Investigation of the axially 
symmetric jet. His results were obtained by employing the "mixing 
length" concept as introduced by Prandtl . His results showed for 
the fully developed jet that 

U c a 1/(Z + e) 

where U is the centerline velocity, Z is the main flow direction, and 
c 

e is the position of the virtual origin measured from the face of the 
nozzle exit along the negative z axis. The above results were based 
on assuming the pressure In the jet was constant. To improve his 
analysis, he considered pressure differences and concluded that the 
pressure differences were so small that no modifications of the veloci- 
ties were necessary. In addition, Tollmein obtained an analytical 
representation of the velocity profile for both the transverse and 
longitudinal component of velocity. The use of one empirical constant 
enabled his analytical results to be brought into very good agreement 
with results of experimental investigations. 

5. Townend, H. C. H., "Flow Induced by a Jet of Air," Rept. No. ARC 1934, 
June, 1934, British Aeronautical Research Council. 

Townend experimentally examined the flow Induced by a jet of air. 

He interpreted his results as consisting of a jet with bounding eddies 
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constituting discrete vortex rings. In addition, It appeared as if the 
jet was alternately accelerated and decelerated as each successive 
vortex ring passed downstream. The overall effect was a kind of pump- 
ing action. 

6. Kuethe, A. M. , "Investigations of the Turbulent Mixing Regions Formed 
by Jets," J. Appl . Mech. , Vol. 11, No. 3, 1935, pp. A87-A95. 

Kuethe considered the turbulent mixing surrounding the jet and 
obtained a solution for the region surrounding the potential core of 
the jet. He assumed no pressure difference between the jet and the 
surrounding air. He used as a first approximation to the longitudinal 
velocity profile in the mixing layer 

U z = (1 - n ,3/2 ) 

where n' = (r - rj)/(r 0 - r^ and r is the radial distance (r.j _<rj r Q ), 
r I is the radius of the potential core, and r Q is the jet radius. 

7. Abramovich, 6. N., The Theory of Turbulent Jets , M.I.T. Press, 
Cambridge, Mass., 1963, pp. 177-208. 

Abramovich presented formulas for the core region and developed 
region of the jet. He employed for the velocity profile in the mixing 
layer surrounding the core the empirical relation 

U z « Uj [1 - (1- n 3/2 ) 2 ] 

where n = (r Q - r)/(r Q - r^ and r Is the radial distance _< r _< r Q ), 
rj is the radius of the potential core, and r Q is the jet radius. This 
relation agreed quite closely with experimental data. Abramovich showed 
that the velocity profiles in the fully developed region were affine. 
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the increase in the thickness or width of the submerged jet (r - r j ) 
was equal to the axial distance (z) times a constant, and the centerline 
velocity of the jet In the fully developed region was inversely pro- 
portional to z. In addition, for the transition he concluded that one 
may assume that the equal velocity lines In the transition region were 
extensions of the equal velocity lines of the potential core region of 
the jet. 

8. Howarth, L., "Concerning the Velocity and Temperature Distributions 
In Plane and Axially Symmetrical Jets," Proc. Cambridge Phil. Soc., 
London , Vol . 34, Pt. 2, 1938, pp. 185-203. 

Howarth made calculations of the velocity profile for the developed 
region of an axisymmetrlc jet. He employed Taylor's mixing length for 
his calculations. Assuming the eddy velocities of the jet were Isotropic 
he obtained results equivalent to employing Prandtl's mixing length. 
Assuming Taylor's mixing length was homogeneous, he obtained results 
equivalent to employing Taylor's modified vorticlty transfer theory. 
Pressure gradients, radial components of velocity, and derivatives of 
radial components of velocity with respect to the axial direction, were 
considered negligible for all calculations. 

9. McElroy, G. E. , "Air Flow at Discharge of Fan-Pipe Lines In Mlnes- 
Part II-Effect of Size and Shape of Pipe and of Adjacent Walls on 
Velocity and Entrainment Ratios," Rept. of Invest. No. 3730, Nov. 
1943, Bureau of Mines, U.S. Dept, of Int. 

McElroy conducted an experimental Investigation and suggested the 
following simple approximate formula for entrainment In the developed 
region of the jet 


E 


m 


I 


C Z/D 


ORIGINAL PAGE IS 
OF POOR QUALITY 


115 


where E m a (Q Q 0 )/Q 0 ** ie ra ^° °f the quantity of entrained air 
In the air stream at any cross section, Q - Q q , to the quantity dis- 
charged Q q and Is called the entrainment ratio, Z is the distance from 
the orifice, 0 is the diameter of the round jet orifice, and C Is a 
constant (about .35). He stated among a number of conclusions: 

(1) The velocity at any point In the air stream and the 
amount of entrained air in the stream at any section vary 
directly with the discharge velocity. 

(2) Conditions that increase centerline velocity ratios 
generally decrease entrainment ratios and vice versa. 

(3) Centerline velocity ratios vary with region. In the 
potential core region they are essentially constant. In the 
transition region, they decrease as the square root of the 
distance; in the developed region, they decrease directly 
with distance. 

(4) The angle of expansion of the developed jet appears 
to be almost 24°. 

(5) Velocity distributions vary with regions. In the 
mixing layer of the core region, the relation of the mean 
velocity to maximum velocity changes rapidly with distance. 

In the transition and developed regions, mean velocities 
probably approximate 40 percent of the maximum velocities. 

(6) Entrainment ratios vary with regions. Analyses of 
the experimental data show continuously increasing entrainment 
ratios with Increasing axial distance. 
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10. Tuve, 6. L. and Prlester, G. D., "The Control of Air Streams In 

Large Spaces," Trans. ASHVE , Vol. 50, 1944, p. 153. 

Tuve and Priester conducted an experimental investigation. They 
concluded that the maximum air velocity at any cross section of the air 
stream beyond 25 diameters downstream varies approximately as follows: 

(1) directly proportional to the exit velocity, 

(2) directly proportional to the diameter of the outlet, and 

(3) Inversely proportional to the distance from the outlet. 

On the basis of an effective outlet area, they obtained an equation 
for entrainment ratio 


£ = 0.785 


R Z /Ae 


/ 


Ae 

wnm 


+ 2 Z Tan 


1.0 


where 

K = constants of proportionality (tabulated) 

R = ratio of maximum velocity to average velocity at any cross 
section (range 2.4 to 3.0) 

e = average spreading angle of the jet (22.5° to 25.4°) 

Z * axial distance from outlet 
Ae = effective outlet area 

11. Squire, H. B., and Trouncer, J., "Round Jets In a General Stream," 
Tech. Rept. No. 1974, Jan. 1944, British Aeron. Research Council 

Squire and Trouncer conducted a mathematical analysis of an axl- 

symmetric jet with a secondary stream. However, within this analysis 

they included the case for a secondary stream of zero velocity. They 

examined the flow in the potential core region and the developed region 

separately and then fitted the solutions together. They determined 
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the inflow Induced by the jet based upon employing cosine distributions 
for the longitudinal velocity profiles. For the mixing layer and the 
developed region, the respective relations were 

s (Uj/2)(1 - cos ir n) [mixing layer] 

U Z a <V 2 W + cos ^ r / r 0 ^ [developed region] 

where U z is the longitudinal velocity, U. Is the exit velocity of the 
jet, U c Is the centerline velocity in the developed region, n = (r Q 

“ r )/(V r I ) and r the ra< ^ a1 position ( rj < r _< r 0 ), rj Is the 
inside radius of the mixing layer, and r Q Is the outside radius of the 
jet. In addition, they made a determination of Inflow velocity by 
assuming a system of sinks along the jet axis. 

12. Cleeves, U. and Boelter, L.M.K., "Isothermal and Non- Isothermal 

Air Jet Investigations," Chem. Engr. Proqr. , Vol. 43. 1947. 
pp. 123-134. 

Cleeves and Boelter conducted an experimental Investigation. They 
concluded that the transverse distribution of the axial component of 
velocity (U z ) In the developed region can be approximately correlated 
by using a dimensionless group. This group was obtained from the 
transfer of momentum In the jet, that Is U z /U c versus r/r M where r^ 

Is the radial position r at U z /U^ * 0.5 and Uq is the centerline velocity 
at the radial position. 

13. Ribner, H. S., "Field of Flow about a Jet and Effect of Jets on 
Stability of Jet Propelled Airplanes," NACA Wartime Report L-213, 
(ARC L6C13, 1946). 

Ribner conducted a mathematical analysis of an axially symmetric 
jet with a secondary flow. He Investigated the Induced flow effects 
on stability and trim of an aircraft as the jet passes near the tall 
surfaces. His analysis was only applicable to the developed jet and 
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he ignored velocity components induced parallel to the jet axis. For 
the spreading of the jet near the origin, he obtained a linear relation- 
ship with the axial distance Z. Far from the origin, he obtained an 
expression for jet spreading as the one-third power of the axial dis- 
tance. In addition, he concluded that the jet-induced flow inclination 
varies nearly inversely as the radial distance from the jet axis within 
the region between the jet boundary and twice the radius of the jet 
boundary at distances greater than 8 orifice diameters downstream of 
the jet exit. 

14. Liepmann, H. W. and Laufer, J., "Investigations on Free Turbulent 
Mixing," NACA TN 1257, Aug. 1947. 

Liepmann and Laufer conducted an analytical investigation of two- 
dimensional free turbulent mixing. They concluded that the mixing 
length theories had lost much of their value in that the main results 
of these theories could be obtained by dimensional reasoning. 

15. Hlnze, J. 0. and Van Der Hegge Zljnen, B. G., "Transfer of Heat and 
Matter in the Turbulent Mixing Zone of an Axially Symmetrical Jet," 
Appl. Scl . Res. , Vol . A1 , 1948. 

Hinze and Van Der Hegge Zljnen made measurements of the radial 
distribution of the mean values of the axial component of velocity, the 
temperature, and gas concentration. They compared these results with 
previous investigators and found that theories based on a constant 
coefficient of shearing stress across the jet flow gave the best agree- 
ment with measured radial velocity distribution in a central zone of 
the jet. 

16. Albertson, M. L., Dal, Y. B., Jensen, R. A., and Rouse, H., 
"Diffusion of Submerged Jets," Proc. of Amer. Soc. of Civil Engjr. , 
Dec. 1948. 
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Albertson, et al . , derived the approximate characteristics of the 
mean flow pattern for flow from slots and orifices. Experimental data 
were used to justify the analytical results and provide the empirical 
constants. In agreement with other Investigators they found it un- 
necessary to make any assumption as to the distribution of turbulence 
to obtain an approximate mean velocity distribution. 

17. Corrsln, S., and Uberol, M. S., "Further Experiments on the Flow 
and Heat Transfer In a Heated Turbulent Air Jet," NACA TN 1865, 
1949. 

Corrsin and Uberoi conducted an experiment on a heated jet issuing 
In still ambient fluid. They suggested that useful results In shear 
flow problems are obtainable with less difficulty by use of Integrated 
equations of motion and reasonable guesses for the velocity profile. 

One conclusion they made was up to a local maximum density ratio, 
p“/p min of about 1.3, simple geometrical similarity still exists In 
the developed region, within the accuracy of measurements. The minimum 
density at a section of the jet was denoted by p min and p°° was the 
density of the receiving medium. In addition, they concluded that 
the pressure and temperature profile functions are basically the 
same as in the constant density jet up to p® / pmin = 1.3. 

18. Corrsin, S., and Kistler, A. L., "Free Stream Boundaries of 
Turbulent Flows," NACA Rept. 1244, 1955. 

Corrsin and Kistler completed an experimental and theoretical study 
of the free stream boundary separating the turbulent fluid from non- 
turbulent fluid. Their suggested theoretical model gave constant mean 
vorticlty at the boundary, local vorticlty production, and a uniform 
suction (induced) velocity. Experimentally they determined for the 
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round jet that the mean velocity at the boundary is chiefly radially 
inward. The angle of spread for the half-cone of the developed region 
was determined to be 10.8° (no-wall). 

19. Phillips, 0. M., "The Irrotatlonal Motion Outside a Free Turbulent 
Boundary," Proc. Cambridge Phil. Soc. , Vol. 51, 1955, pp. 220-229. 

Phillips analytically considered the irrotational motion of an 
infinite fluid when the normal velocity across a plane is a stationary 
random function of position. He postulated conditions which corres- 
ponded closely to the motion outside a free turbulent boundary. 

20. Stewart, R. W., "Irrotatlonal Motion Associated with Free Turbulent 
Flows," J. Fluid Mech. , Vol. 1, 1956, pp. 593-606. 

Stewart conducted a theoretical examination of the irrotatlonal 
motion external to a self preserving turbulent wake and jet. He found 
a mean flow towards the center of the jet and pointed out that it is 
possible for the vorticity free fluid between bulges of turbulent fluid 
to obtain the mean velocity of the turbulent fluid. However, he con- 
cluded that this was a very short range effect and the motion in the 
external region can be considered Irrotational. From similarity consid- 
erations, Stewart expressed the mean flow near the axi symmetric jet 
boundary as U = Cp“\ where p are the colatitude angle and the 
radius in a spherical coordinate system with the jet exit located 
at the origin and directed along the z axis. For a jet emerging nor- 
mally through an infinite plane baffle plate, he found that the velocity 
potential <t> can be expressed by 


♦ ■ C tan zn (p sin $) 
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where = the half cone angle containing all of the turbulent region 
of the jet (about 12.5°). 

21. Miller, D. R. and Comings, E. W. , "Static Pressure Distribution In 
the Free Turbulent Jet," J. Fluid Mech .. Vol. 3, 1957, pp. 1-16. 

Miller and Comings made measurements of mean velocity, turbulent 

stress, and static pressure in the mixing layer of a jet of air issuing 

from a slot (2-D) Into still air. Appreciable deviations from constant 

pressure conditions were found. They obtained negative static pressure 

readings everywhere in the mixing layer except In the potential core 

wedge. 

22. Ricou, F. P. and Spalding, D. B. , "Measurements of Entrainment by 
Axl symmetrical Turbulent Jets," J. Fluid Mech., Vol. 11. 1961. 
pp. 21-32. 

Ricou and Spalding made measurements of entrainment by means of 
a porous-walled, cylindrical chamber. The flow rate through the porous 
v/all was adjusted until no axial pressure gradients could be detected. 

The entrained fluid as a result of the radial Inflow was then presumed 
to be equal to that of a free jet. For the fully developed jet they 
obtained a relationship for the mass entrainment as a function of axial 
position. As In previous work by other Investigators, they had to 
determine a numerical constant (Kj of Equation (35)). Their experimental 
value was 0.282 which was within the range of values 0.22 to 0.404, 
determined by earlier Investigators. 

23. Uygnanski, I., "The flow Induced by Two-Dimensional and Axlsymmetrlc 
Turbulent Jets Issuing Normally from an Infinite Plane Surface," 
Aeron. Quart ., Vol. XV, 1964, pp. 373-380. 

Wygnanskl made an analytical investigation of turbulent jets. 

The pressure distribution on the surface from which the jet was issuing 
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was determined. In his analysis, he replaced the jet with a line of 
sinks with variable strength. His results agreed with the experimental 
data he presented. 

24. Gentry, C. L. and Margason, R. J., "Jet-Induced Lift Losses on VTOL 
Configurations Hovering In and Out of Ground Effect," NASA TN 
D-3166, Feb. 1966. 

Gentry and Margason conducted an experimental Investigation of 
pressure distribution on a plate from which a jet flow was exhausting. 

The case where the plate was flush with a large wall and the case 
where the plate was exposed to the atmosphere on both sides were 
studied. Pressure coefficients on the plate were measured In both 
cases. In comparing their results for the large wall with those of 
Wygnanskl, they show that Wygnanskl's data Is about 25 percent less. 

In general, many mathematical analysis and experimental studies 
of jets issuing from nozzles with diameters ranging from 0.102 to 20 
Inches and with velocities ranging from 13 to 1000 ft/sec have been 
made. For the developed region, using the concept of similarity, 
theories for mean velocity distributions give about the same picture. 
However, they yield a velocity distribution which Is too sharp an apex 
on the jet axis In comparison with experimental results. As Is to be 
expected, near the jet boundary comparison of the analytical velocity 
profile with the experimentally determined distribution show a devia- 
tion from the measured distribution. 
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APPENDIX B 


Derivatio n of the Equations for the Aerodynamic 
Field of an Axl symmetric Jet 

This Appendix contains derivations of the equations for determining 
the aerodynamic velocity field for an axlsymmetrlc jet according to the 
method which make use of the Stokes potentials. The derivations are for 
a semi -Infinite domain where Green's function Is given by 2/R. Starting 
from Equation (15) 

U = U ST + E p 05) 

where 

Ust = * V 1 x U l> x “J dT l 

and 



V(|) u 


1 


dS, 


Since the flow is axisymmetrlc (v^ x 1 j), 

^ 0I " ^ST C3n Wr ^^ en aS 




^gj X R 


61 

T 


dr. 


the equation 
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For the geometry of Fig. 11, noting that e may be put equal to zero 
without loss of generality, R can be written as 


R » (r cos 0 ] - r^ 1 fl - r sin G 1 \ Q} + (z-z^ i 2l 

The source point coordinate system Is the appropriate system for the 
volume integration. Substituting this result into the equation for 

u st 


St ” 


kf T 


( z * z l) Vl " n el ( r cos e i" r i) 1 


l'l' 'zl 


dt. 


which can be separated Into r and z components In the field point 
coordinate system. 

2it °° r. 


U 


Sxr “ 4? 


M S I 


1 cos 0.| fi 0 i (z-z-j ) ri drj dZi de^ 
R 3 " 


0 -« o 


Stz 


2 tt °° 

Hi 

w o --•'0 


1 ^ 01 (r^r cos 0^ 


ri dri dZi dO^ 


Making use of the Identity that 




? 


the equation for E p may be written as 


E P * h ( ^j ) u i * dS i 
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Using the geometry of Fig. 12, noting that 6 may again be put equal to 
zero, R can be written as 


R = (r cos e 1 - r ] ) 1 fl - r sin 0, i 0l + z i 


zl 


Substituting this result into equation for F along with U. for 

P J 

U-j • n yields 


; . u j f (r cos 9 l - r i> Vl - r s1 " e i <ei + 2 ’ 2 i ... 

» 2, J 7 — dS l 


Separating the above relation into r and z components yields 


E r — i 

pr 2 tt 


ff 


U.. f 2lT r r l (r-r^ cos 0 1 ) 

—3 r l dr l de l 


^if 2 V ri r, 

pz = 2T 3 dr l de l 

J o Jo K 


Consider next the Integral over 6 1 for U STr 


\ 


2ir 


cos 6-| dO^ 
0 (a-b cos e 1 ) 


m 


where a = r 2 + r 2 + (z-z^ 2 and b = 2^ 


Let 


t, -2 


j: 


de. 


(a-b cos e-j) 


3/7 
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then 


I 




cos e-j de-j 
(a-b cos e.j)^ 2 


a I^b I 


2 



de, 


(a-b cos e-j) 


V? 


Solving this relation for I 2 , the integral over e-j for U^, yields 


I = — 
*2 b 


J 


de. 


q (a-b cos e^) 


3/T 


2 

b 


1 


de. 


0 (a-b cos e-j) 


7m 


These Integrals may be evaluated in terms of complete elliptic integral 
24 

functions. Thus, 


\ 


de. 


0 (a-b cos e^) 


377 


A E < A ’ k > 


where 

g 3 2/ /a + b, k 2 3 2b/ (a + b) 
and 

A 3 sin" 1 / b (1-cos o 1 )/k 2 (a-b cos 0-j) 

Substituting the value for e 1 In the relation for A yields A = sin 
1 3 n/2. Thus, E (A,k) 3 E U/2, k) 3 E (k) a complete elliptic 
Integral of the second kind. Now, substituting for g, one obtains 


2a 

EP 




(a-b cos e<j) 


. 4a 

377" FTa=FT 


UkL 

/a + b 
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Similarly, 


2 

5" 


I 


de 


0 (a-b cos e-j) 


’ . 4 . K (k)„ 

b / a + b 


where K (k) Is a complete elliptic Integral of the first kind. Sub- 
stituting these results for the right side of the equation for I 2 yields 

Fb E (k) " K (k) 


*2 = 

b / a + b 


Likewise Integrating over e-j for U<. with 




J 


do. 


g (a-b cos o-j ) 


77? 


r - 2 


I 


cos o.| de-j 


g (a-b cos o-j) 


37? 


yields 


Io = 


2 /a + b (a-b) 


4 e (k) - — 4 = [A-e (k) - k (k)] 


b 4Tb la ' b 


Substituting the respective results Into the Integrand of the equations 
for U $Tr and U $TZ , yields 


J = — 
J Sxr 2irr 


f n e, (z - z l> 

" S <7 (r + r, 12 + 


S /(r + rj)2 ♦ (z-zj)2 


^ r 2 + r 2 + ( z — z i ) 2 
(r-^) 2 + (z-z^ 2 


E 00 - K (k) 


dr-| dz-| 


Stz 


L. f 

2 ’J S 477777 


S /(r + r,) + (z-z,)' 


r 2 -r\ + (z-z,) 2 


*1 


(r-r-j ) 2 + (z-z^ 2 


(40) 
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E (k) - K (k) 


drj dzj 


(41) 


In order to obtain Equations (42) and (43), the equations for E 


Pr 


and E pz are rewritten as 


E * 
pr 


U.r f 2 * f r l r 1 dr-j de ] U, f £1T C'l r^ cos e 1 

** Jo Jo (a " b cos 0 i )3/2 " Jo Jo (a ' b cos e i )3/2 


2tt *r, J. 


dr-j de-| 


'P2 


. V f 2 *f ri r i dr i de i 

2 " Jo Jo ( a " b cos e i) 


3/2 


Integrating over e-j, yields 


r r 

2 2 2 ) 

r -ri-z 

l ! 

E (k) + K (k) 

*r /(r + r-j) 2 + z 2 

L (r-r,) 2 + z Z 


dr. 


(42) 


: = 2 U j 1 f ^ 

PZ Jr /(r + r, ) 2 + 


1_M 


2 . .2, dr l 


(43) 


/(r + r-j ) 2 + z 2 + z ] 
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APPENDIX C 


Algorithms for Complete El 1 i pti c Functions 

The following algorithm computes the complete elliptic integral of 

12 

the first kind with modulus k . 

t 

K(k) = {[(0.032024666t + 0.054544409)t + 0.097932891]/ 

1.3862944} - {[ (0 . 0 109449 1 2 1 + 0.060118519 )t + 


0. 12475074] t + 0.5 } ln(t) 

2 

where t = 1 - k 

The following algorithm computes the complete elliptic integral of 

1 3 

the second kind with modulus k . 


E(k) = {[0.O4O9O5O94t + 0.085099l93)t + 0.44479204]t + 

1.0} - {[(0.01 382999t + 0.08150224) t + 0.24969795]t} ln(t) 

2 

where t = 1 - k 



